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MODELING  AND  CONTROL  OF  ALGAE  HARVESTING,  DEWATERING 
AND  DRYING  (HDD)  SYSTEMS 


Abstract 

By 

FENGMING  LI 

Alternative  energy  is  a  very  important  topic  since  there  are  some  issues  with  the 
continued  use  of  fossil  fuel  as  a  major  source  of  energy.  The  alternative  energy 
studied  in  this  research  is  algae,  which  has  received  a  lot  of  attention  because  it  is 
easy  to  cultivate,  the  harvest  cycle  is  only  around  10  days,  the  yield  is  5000~10000 
gallons  per  year  per  acre,  the  oil  extracted  can  be  used  directly  in  a  variety  of 
applications,  and  there  is  no  further  investment  required  in  building  transportation 
facilities.  The  key  barrier  to  making  algae  oil  a  commercial  success  in  the  fuel  market 
is  the  cost  of  production.  One  of  the  major  costs  of  extracting  oil  from  algae  is  the 
drying  process,  which  is  about  30~50  %  of  the  total  processing  cost.  In  this  research, 
experiments  are  conducted  and  a  predictive  model  that  can  be  used  to  develop 
optimal  operating  strategies  for  HDD  systems,  which  provide  the  most  efficient 
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method  for  dewatering  algae.  The  predictive  model  is  developed  using  fluid 
dynamics  principles,  and  the  dynamic  aspects  of  the  web  transportation  system. 
Simulation  results  show  that  the  predictive  model  can  be  effectively  used  to  operate 
the  HDD  system  in  different  situations. 


XVI 1 1 


Chapter  1  Introduction 


1.1  Energy  Crisis  and  Climate  Change 

A  major  turning  point  occurred  in  the  mid-1700s  during  the  industrial 
revolution.  The  industrial  revolution  did  not  only  change  the  way  people  live, 
but  also  had  an  effect  on  the  Earth's  ecology.  Fossil  fuels  became  a  major  source 
for  power,  and  it  has  been  widely  used  since  then.  There  is  no  doubt  that  fossil 
fuels  played  an  important  role  in  industrialization,  but  they  also  had  some 
adverse  effects  on  the  environment.  According  to  lEA's  report  [1],  around  81.4% 
of  the  world's  energy  was  generated  from  fossil  oil  in  2007.  In  other  words,  fossil 
fuels  dominate  the  energy  supply.  With  the  rising  population  and  increasing 
demand  of  energy,  fossil  fuels  will  be  depleted  since  it  is  a  finite,  nonrenewable, 

resource. 

Another  issue  created  by  fossil  fuel  is  the  greenhouse  effect,  where  the 
balance  between  radiation  entering  and  leaving  the  earth  has  been  disrupted. 
Figure  1.1  shows  a  representation  of  the  exchange  of  energy  between  the  Sun, 
the  Earth's  surface,  the  Earth's  atmosphere,  and  outer  space  [2].  Extra 
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greenhouse  gas  has  been  introduced  into  the  system  over  the  years  since  fossil 


fuels  were  used  broadly.  It  leads  to  an  increase  of  radiation  absorbed  by 
greenhouse  gas,  and  the  original  balance  is  disturbed.  As  a  matter  of  fact, 
according  to  the  report  [3]  from  the  Intergovernmental  Panel  on  Climate  Change 
(IPCC),  the  levels  of  Carbon  Dioxide,  Methane,  and  Nitrous  Oxide  have  risen 
considerably  since  the  Industrial  Revolution.  Another  report  [4]  issued  by  the 
IPCC  demonstrated  that  the  increase  of  global  average  temperature,  the 
increase  of  global  average  sea  level,  and  the  decrease  of  northern  hemisphere 
snow  cover  may  be  related  to  the  greenhouse  effect.  The  report  also  shows  the 
upward  tendency  of  the  global  average  surface  temperature  and  global  average 
sea  level  seem  to  match  the  upward  swing  of  greenhouse  gases.  If  that  is  true,  it 
will  be  difficult  to  reverse  the  trend  of  increasing  global  average  surface 
temperature  and  sea  level  unless  we  reduce  the  use  fossil  fuels.  Besides,  report 
provided  by  IPCC  [4]  illustrates  the  impact  of  the  global  average  temperature 
change  in  the  period  1980-1999.  It  shows  the  impacts  are  all-round  us,  and  the 
result  will  be  fatal  if  we  don't  respond  quickly. 
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Fig.  1.1  Simplification  of  the  flows  between  space,  the  atmosphere,  and  the 

Earth's  surface  [2] 


1.2  Alternative  Energy 


Because  of  the  energy  crisis,  and  the  fatal  side  effects  of  fossil  energy, 
searching  for  reliable  and  never-depleting  alternative  energy  is  an  important 
and  urgent  issue  nowadays.  There  are  some  common  types  of  alternative 
energy,  such  as  solar,  wave,  biomass,  tidal,  and  wind  power,  discussed  in  the 
past  decades.  Among  them,  bio-fuel  is  one  of  the  most  popular  alternative 
energy  resources  for  the  following  reasons.  First,  major  modifications  are  not 
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necessary  to  switch  to  the  use  of  bio-fuel.  Second,  it  can  reduce  greenhouse  gas 


emissions,  and  sustain  longer  than  traditional  fuel.  Third,  it  is  a  renewable 
source.  Fourth,  bio-fuels  are  easily  available  from  common  biomass  sources, 
environmentally  friendly,  and  biodegradable  [5][6]. 

Among  all  kind  of  biomasses,  algae  is  one  of  the  most  potential  species  and  it  is 
arising  scientific  and  public  attention  due  to  the  reason  that  it  is  low-input,  high-yield, 
up  to  30  times  more  energy  per  acre  than  land  crop  be  produced  is  claimed  based  on 
the  experiment,  to  produce  bio-fuels,  it  allows  solar  energy  to  be  stored,  and  it  can 
be  used  in  today's  engine  and  transport  infrastructure  just  like  other  bio-fuel. 
Moreover,  unlike  other  renewable  energy  resources  such  as  wind,  tidal,  and  solar 
power,  algae  can  be  cultivated  and  harvested  systematically.  After  dewatering,  oil  in 
algae  can  be  extracted  to  produce  bio-diesel,  and  the  residue  can  be  used  as 
additives  in  animal  feed  or  used  as  feedstock  in  fast-pyrolysis  process  to  produce 
bio-oil.  Besides,  algae  are  capable  of  making  CO2  as  a  carbon  source  just  like  other 
plants.  Figure  1.2  shows  the  process  of  closed  carbon  cycle  of  algae.  Besides,  they  are 
one  of  the  most  photosynthetically  efficient  plants  on  the  earth,  and  they  use 
photosynthesis  to  convert  solar  energy  to  chemical  energy.  The  energy  is  storied  in 
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the  form  of  oils,  proteins  and  the  others.  The  algae  oils  can  be  made  into  biodiesel, 


and  biodiesel  is  gaining  public  and  scientific  attention  since  it  is  competitive  with 
traditional  petroleum  diesel  fuel  technically  [7],  Algae  also  have  some  other 
advantages.  They  can  be  cultivated  in  either  marine  or  fresh  waters  so  that  the 
occupancy  of  land  can  be  avoided.  And  as  you  may  know,  over  70%  of  the  surface  of 
the  earth  is  water.  In  other  words,  algae  have  plenty  of  places  to  grow.  Moreover,  the 
harvesting  cycle  of  algae  is  about  1-10  days.  It  is  way  beyond  fast  to  collect  bio-oil.  In 
short,  algae  have  a  lot  of  potential  to  be  substitute  for  fossil  fuel.  Therefore,  algae 
can  be  expected  to  be  the  major  fuel  source  in  the  future. 


Fig.  1.2  Closed  carbon  cycle  of  algae 


1.3  Algae  oil  extraction 


There  are  a  lot  of  reasons  to  develop  algae  oil  technology.  However,  some 


technique  issues  still  need  to  be  overcome.  One  of  the  issues  is  that  algae  oil 
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extraction  is  a  very  costly  process.  Hence,  further  research  on  how  to  efficiently 


extract  oils  from  algae  is  necessary.  In  this  study,  a  new  algae  harvesting  system 
will  be  proposed,  with  a  focus  on  mathematically  modeling  and  control  the 
system. 

1.3.1  Proposed  HDD  System 

ARPA-E  provided  funding  to  a  team  lead  by  Univenture  /  Algaeventure 
Systems  (AVS)  to  develop  an  innovative  technology  for  removing  the  water  from 
suspended  algae  with  dramatically  reduced  energy  consumption  by  utilizing 
surface  physics  and  capillary  action  to  more  effectively  harvest,  dewater,  and 
dry  algae  (HDD).  Highly  efficient  HDD  offers  the  potential  to  transform  the 
economics  of  algae-based  bio-fuel  production,  removing  the  major  barrier  to 
large-scale  commercialization  of  this  renewable  alternative  fuel  source.  A 
lab-scale  prototype  system  at  AVS  has  demonstrated  95%  energy  reduction  in 
comparison  to  conventional  centrifuge  methods.  Based  on  surface  physics  of 
materials,  it  uses  several  interacting  web  surfaces  to  continually  extract  algae 
from  very  low  algae-water  concentrations  and  dry  it  to  less  than  5%  moisture 

content. 
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Harvesting  and  dewatering  algae  has  long  been  cited  as  the  major  hurdle 


to  making  algae  a  commercially  viable  feedstock  for  bio-fuels  production. 
Industry  expert  John  Benemann  noted  the  logistical  difficulties  associated  with 
algae  harvesting  and  dewatering  as  early  as  1982,  and  recently  (2008) 
re-asserted  that  "at  present,  there  are  no  low-cost  harvesting  technologies 
available".  This  remains  a  serious  deficiency  in  commercial  systems  today. 
Dewatering  processes  are  important,  because  as  the  sludge  becomes  drier,  the 
lipid  extraction  process  becomes  more  efficient.  Transforming  raw  solution  to  a 
dense  sludge  sufficient  for  lipid  extraction  has  historically  been  a  multistage 
energy  inefficient  process  consuming  30%  -  50%  of  the  total  cost  of  algae 
cultivation. 

1.3.2  HDD  Process 

The  proposed  HDD  (harvesting,  dewatering,  and  drying)  process  takes 
advantage  of  water  surface  characteristics,  wherein  water  adhesion  and 
capillary  action  is  used  to  efficiently  separate  the  sparse  algae  colonies  into  a 
dry  cake-like  biomass  with  extremely  low  moisture  content  that  can  be  directly 
processed  for  oil  extraction.  The  proposed  HDD  Laboratory  System  (Figure  1.3) 
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consists  of  an  input  load  well  that  accepts  a  steady  flow  of  water  containing 


low-to-moderate  concentrations  (0.3  to  40  g/L)  of  algae  biomass.  A  thin  film  of 
algae  is  collected  on  a  porous  filter  web  that  passes  under  the  input  load  well. 
Gravity  causes  a  portion  of  the  water  to  drain  through  the  porous  filter  while 
the  filter  pore  size  is  small  enough  to  trap  algae  and  transport  it  to  contact  the 
capillary  action  web  traveling  in  the  opposite  direction.  The  parallel  moving 
capillary  belt  will  cause  water  to  wick  away  from  the  algae  using  capillary  action 
where  water  contacts  the  hydrophilic  surface.  This  unique  capillary  action 
provides  excellent  water  removal  without  differential  pressure  that  could 
embed  the  algae  into  the  filter  and  lead  to  filter  pore  fouling.  A  final  drying 
stage  is  typically  an  evaporative  approach,  and  the  lab  model  consists  of  an 
electrically  heated  plate  that  is  in  close  contact  with  the  filter  membrane  to 
remove  residual  surface  moisture.  Water  draining  from  the  system  may  be 
reused  as  culture  media  or  used  for  other  downstream  processing  steps. 
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Separates  algae  from  water,  dewaters,  dries 
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Fig.  1.3  Proposed  HDD  system  schematic 


1.3.3  Laboratory  System  Performance 

Compared  with  other  dewatering  systems,  the  HDD  device  does  not  put 
energy  into  the  system  such  as  to  maintain  differential  pressure,  operate  a 
centrifuge  or  bubble  air  through  the  solution  for  subsequent  skimming,  or 
pressure/flow  based  filtering.  The  only  energy  needed  is  to  operate  a  low-drag 
transport  belt  and  a  final  stage  heater,  providing  transformational  production 
cost  savings  when  compared  with  conventional  methods  (i.e.  centrifuge). 
System  performance  has  been  tested  under  controlled  conditions  for  a 
representative  set  of  Algae  (Table  1.1). 
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Process  runs  were  performed  over  a  fixed  processing  time  (10  minutes), 


belt  rate,  and  heating  under  the  following  constraints  and  limitations: 

•  System  performance  requires  a  thin  layer  of  solids  (the  filter  cake)  to  be 
"painted"  onto  the  filter  belt  to  facilitate  filtration  and  subsequent  drying. 

•  Raw  effluent  throughput  is  limited  to  a  flow  rate  that  maintains  filter  cake 
consistency. 

•  Suboptimal  "Off  the  Shelf"  filters  were  tested. 

Results  showed  that  the  energy  needed  to  dewater  to  20%  solids  was  6.7 
Whrs  in  10  minutes  based  on  motor  current,  while  the  energy  needed  to 
dewater  and  dry  to  <5%  moisture  was  206  Whrs  in  10  minutes  based  on 
motor/heater  current.  Furthermore,  the  AVS  HDD  Lab  System  running  at  500 
L/hr  for  a  3  g/L  algae  concentration  used  only  40  watts-hrs  for  dewatering  1498 
g  algae  (dry  weight).  The  resulting  energy  consumption,  0.0267  Whr/g,  is  94.7% 
less  than  the  best  post  flocculation  centrifugal  dewaterer  (0.5  Whr/g)  while 
doubling  the  solids  concentration  (20%  for  the  proposed  HDD  vs.  10%  for  the 
centrifuge).  Beyond  the  energy  cost  savings  for  the  dewatering  stage  of 
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processing,  this  methodology  also  provides  the  opportunity  for  improvements 


in  overall  process  performance: 

•  Less  drying  (Figure  1.4):  The  doubling  of  the  algae  solids  concentration 
during  dewatering  reduces  the  amount  of  energy  needed  by  the  heating 
stage.  Indeed,  when  optimized,  we  believe  the  HDD  can  achieve  25%  solids 
by  weight  from  the  dewatering  stage  alone.  Based  on  a  latent  heat  of 
vaporization  model  and  conservatively  assuming  Lwater=  800  kcal/kg,  3.5 
Wh/g  will  be  needed  by  the  AVS  HDD  to  drive  75%  water  (250  g/L  algae 
concentration)  down  to  5%  water  (950  g/L  algae),  while  the  centrifuge 
system  will  need  8.32  Wh/g  starting  at  80%  water  (200  g/L  algae).  This 
represents  a  56%  savings  in  the  drying  stage!  Combining  both  dewatering 
and  drying,  the  energy  savings  is  approximately  58%  (3.7  Wh/g  AVS  HDD  vs. 
8.82  Centrifuge).  Because  harvesting,  dewatering,  and  drying  are  estimated 
to  consume  30%  to  50%  of  total  algae  oil  production  costs,  the  proposed 
HDD  may  represent  a  savings  of  between  17%  and  29%  of  total  algae  oil 
production  costs.  However,  these  numbers  represent  theoretical  savings. 

•  Facilitates  drying:  Thin  "painted  on"  belt  algae  maximizes  the 
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surface  area  to  volume  ratio  facilitating  downstream  drying. 


including  the  possibility  of  minimizing/omitting  the  heating 
element. 

•  Minimizes  cell  analysis:  Pressurized  belt  filtration  system  analyses 
algal  cells  releasing  valuable  oil,  which  will  escape  with  the  water.  A 
low  pressure  HDD  design  would  eliminate  this  issue. 

•  Increased  filter  belt  life  expectancy:  Beyond  the  benefits  of 
minimized  fouling,  vendor  feedback  indicates  low-pressure  design 
extends  belt  life  due  to  reduced  stretching. 


12 


Table  1.1:  AVS  HDD  Laboratory  System  Performance.  Data  is  presented  as 

average  across  all  runs. 


Algae 

Botryococcus 

braunii 

Chlorella 

vulgaris 

Euglena 

gracilis 

Nannochlorops 

is  salina 

Pond  Water 

Marysville,  OH 

Image 

M 

1  ^7^ 

Size  Range  (i-im) 

15-35 

2-10 

15-40 

0.5-2 

0.5-40 

Density  (g/L) 

2 

2.8 

3.1 

40.2 

1 

Algae  Energy  Content 

(Watt-hrs/gram  Algae) 

7.2 

5.2 

5.6 

6.8 

4.6 

#  of  Process  Runs 

2 

IS 

12 

39 

53 

Flow  rate  (L/hr) 

456 

342 

342 

60 

510 

Escaped  Algae  (%)* 

0.5 

8 

1 

35 

3 

Residual  Water  in  Dry 

Flake  (%) 

1.2 

0.5 

1.4 

0.4 

1.0 

Dry  Solids  Yield  (g) 

152 

147 

175 

261.3 

83 

Escaped  Algae  are  the  percentage  of  Algae  that  transferred  with  the  water  through  the  filter. 
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Fig.  1.4  Dryer  stage  -  Energy  consumed  when  driving  starting  concentration 
to  5%  water  based  on  latent  heat  of  vaporization 
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Fig.  1.5  Micrographs  of  dry  algae  filter  adherence  post  AVS  HDD  (left)  and 
Pressure  Filtration  (right)  Magnification:  Top-left  170x;  -right  lOOx;  Bottom-left 

650x;  -right  850x 


1.3.4  Objectives 


A  prototype  of  the  proposed  HDD  system  was  constructed.  Dynamic 
models  of  the  HDD  process  will  be  developed  to  capture  the  fundamental 
physics  that  are  core  to  the  operation  of  HDD.  Physical  models  of  core 
phenomenon  such  as  filtering,  capillary  action,  drying,  and  water  cohesion  (e.g. 
surface  tension)  will  be  established  with  suitable  fidelity  to  permit  combining 
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these  models  into  a  cohesive  integrated  functional  model.  Analytic  models  will 


be  developed  to  enhance  ease  of  use  and  form  the  foundation  for  an  integrated 
approach  to  system  and  control  design. 

The  laboratory  HDD  system  will  be  used  to  develop  both  a  parametric 
model  of  equipment  operation  and  a  first  order  analytic  model  of  the  drying 
process.  The  two  models  will  be  combined  to  predict  the  moisture  content  of 
the  dried  algae  output.  The  thesis  is  organized  as  follows:  introduction  in 
chapter  1,  mathematical  formulation  and  CFD  method  in  chapter  2,  web 
modeling  and  tension  control  in  chapter  3,  experimental  methods  in  chapter  4, 
results  and  discussion  in  chapter  5,  and  future  work  in  chapter  6. 

1.3.5  Contributions 

In  this  research,  a  predictive  model  that  can  be  used  to  operate  the  HDD 
system  is  developed.  Using  the  model,  operators  can  predict  the  fluid  dynamics 
in  the  system,  and  decide  on  the  "optimal"  operating  settings  for  both  the  web 
transport  system,  and  the  heater  to  reduce  the  overall  operating  cost.  Moreover, 
the  methods  used  in  this  research  can  also  be  adopted  to  simulate  the  scaled-up 
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system  and  can  be  applied  to  different  geometries,  filters,  and  different  kinds  of 
algae.  Hence,  the  cost  of  trial-and-error  experiments  and  the  development  time 
for  the  system  can  be  reduced.  Overall,  not  only  does  this  research  have  the 
potential  to  reduce  operating  costs  save  and  time  but  also  offers  a  new 
approach  to  the  study  of  HDD  systems. 
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Chapter  2  Mathematical  Modeling  of  the  HDD 


According  to  observation  and  experimental  data,  some  important  physical 
processes  involved  in  the  HDD  are  studied.  The  objective  of  this  chapter  is  to 
establish  the  governing  equations  and  create  a  simulation  model  for  each  part  of 
HDD  system  based  on  the  physical  laws. 

2.1  Transport  equation 

Generally,  the  scalar  transport  equation  can  be  stated  as  follows: 

fc.v.  ^  ■  ndA  =  ^  r^VTi  ■  ndA  +  ^  S^dV 

(2-1) 

Where  q  is  the  scalar  quality,  V  is  volume,  c.v.  is  a  control  volume,  p  is  density, 
n  is  the  outward  point  normal  vector  on  the  control  surface,  A  is  the  surface  area, 
and  V  is  the  gradient  operator,  priv-fi  is  the  convective  vector  normal  to  the 
surface,  and  F^Vr)  ■  fi  is  the  diffusive  vector  normal  to  the  surface.  The  net  flux 
passing  though  the  boundary  of  control  volume  can  be  expressed  in  the  following 
form: 


17 
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(1-1) 


Where  fi  is  diffusive  or  convective  vector  normal  to  ith  face,  and  k\  is  the 


surface  area  of  the  ith  face. 


2.2  Discretization  method 


2.2.1  Discretization  of  diffusive  term 


For  computational  convenience,  the  transport  equation  is  discretized.  If  the 


geometry  is  two  dimensional,  from  equations  (2-1)  and  (2-2),  the  discretization  of 


the  diffusive  and  source  terms  can  be  stated  in  the  following: 


f  r.vii  •  fidA  +  s,dv  =  {r,A  1)^  -  {r,A  +  (r,A  - 


+  SrjAV 


(2-3) 


Where  the  gradient  of  rj  can  be  approximated  using  linear  profile 


approximation. 
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Because  the  average  value  of  S,,  may  be  a  function  of  the  dependent  variables, 
it  can  be  approximated  in  the  following  form: 


Srj  ~  +  Sjj  pTjp 


(2-4) 


Hence,  the  discretization  of  the  diffusive  and  source  terms  can  be  formulated 
in  the  following: 


r  A 

^rj.e^e 


Ve-Vp 

Xp-Xp 


r  A 

^rj.w^w 


Vp-Vw 

Xp—Xy\r 


"h  {,^ri,u  ^ri,pV 


(2-5) 
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2.2.2  Discretization  of  the  convective  terms 


From  the  discussion  above,  the  values  of  scalar  variables  are  located  at  the 
centers  of  each  cell.  However,  the  values  on  surface  are  required  for  the 
convection  term.  In  this  study,  face  values  of  scalar  variable  are  function  of 
values  at  the  center  of  the  cells.  Here,  Quadratic  Upwind  Interpolation  for 
Convective  Kinematics  (QUICK)  [14]  and  second-order  Upwind  [15]  schemes 
are  adopted  for  this  study.  There  are  three  cell  points  involved  in  the  scheme.  If 
the  control  volume  is  as  shown  in  Figure  2.1,  the  discretized  form  can  be 
expressed  as  follows: 


Ue  = 


A  [  r|p  H — +  (1  ~  A)  [■ 

[dE+dp  dE+dpx  ^  ^  [ 


dw+2dpx 


dfi+dpx 
dfi 


A  +  T— ^Uee  +(1-A)  ■ 

LdEE  +  dE  dEE  +  dE 


dw+dp 
dpx+Zdp 
dpx+dp 


Tip 


^Px 


Ue 


dfi+dpx 

dfi 


dEE+dp 


]  if(v  ■  n)e 
if(v  ■  n), 


>  0 


<  0 


(2-6) 
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Py 
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dpy  -l-  2df, 
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Py 
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-tin 


dpy  -l-  d] 
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(2-7) 
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Fig.  2.1  Two  dimension  cell  structure 

The  value  of  X  is  setting  to  6/8  for  the  QUICK  scheme.  However,  switching 
between  schemes  is  made  possible  by  changing  X.  For  example,  if  X  is  1,  equation  5 
will  be  a  second-order  central  interpolation.  If  A.  is  0,  a  second-order  upwind 
scheme  will  be  developed. 

With  the  approximation  of  scalar  variables  on  the  surface,  the  convective 
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terms  can  be  obtained  as  follows: 


ipv^r]A)e  -  (pVxVA)^  +  (pVy-riA)^  -  {pVyViA)^ 


(2-8) 


Where  rj  can  be  expressed  by  using  equations  (3-7)  and  (3-8). 


2.2.3  Discretization  of  transient  terms 


In  the  transient  problem,  time  is  another  "dimension"  needs  to  be  considered. 
However,  discretization  in  the  temporal  is  different  than  in  the  spatial  dimensions. 
The  main  difference  is  that  the  time  dependent  variables  only  influence  the  system 
in  the  future.  In  other  words,  there  is  no  backward  affect  in  this  domain. 
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(2-9) 


From  the  previous  discussion,  the  spatial,  diffusive,  convective  and  source 
terms,  have  an  appropriate  discretization,  thus  the  equation  above  can  be 
simplified  into  the  following  form: 

Or: 

^  =  f(Kt)) 

Here,  the  temporal  discretization  can  be  expressed  as: 

First  order:  Vm+i-Vm  ^  p.  . 

At  ^ 

Second  order:  ^ 

Where  i/m-i  is  the  value  of  the  scalar  at  the  previous  time  step,  i/m  is  the 
value  at  current  time  step,  and  i/m+i  is  the  value  at  the  next  time  step. 
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2.2A  Discretized  equation  of  the  system 


If  the  scheme  is  applied  to  the  system  of  equations,  scalar  variables  at  the 
surface  can  be  expressed  in  the  following  form: 


^e,m 


kiTlE,m  +  ^2^P,m  +  k3Tlw,m  if(v  '  n)e,m  >  0 
k4TlEE,m  “I”  ks^lE,!!!  k5T|p  m  if(v  ■  n)e  <  0 


(2-10) 


Where 


ki  = 

k4  = 


A.dp 

dp+dp 

dEE+dfi^ 


ko  = 


kr  = 


dE+dp 

dE+dEE 


+  (1-A) 


dvv+2dp 

dw+dp 


k3  =  (X-l) 


+  0-A)^,  k,  =  (X-l) 


d^-l-dp 

dE 

dEE+dE 


lln,ni 


k7riN,m  +  kgTlp^m  +  k9r|s,m  if(v  ■  n)n  ,n  >  0 
kiollNN,m  +  kiir|N,m  +  kizTlp^m  if(v  '  n)n  m  <  0 


(2-11) 


where 


k? 


A.dp 

dfj+dp 


ks 


^dN 

djsj+dp 


+  (1-A) 


dc+2dp  I 


(A-l) 


dp 

dfj+dp 
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ki  n  — 


Adfj  ,  dp+2dN 


"10  -  j  ^11  -  H  +rl 

Onn  +  Un  Ufj  +  UiviN 


+  (1-A) 


dp+dfj 


ki2  —  1) 


dN 


dNN+^N 


If  we  combine  all  discretization  equations  together  into  equation  (2-1),  the 
following  form  can  be  obtained: 


^ P,mVp,m  ^W,mllw,m  ^E,mllE,m  ^N,mllN,m  ”1"  ^S,mlls, 


m 


-I-  S'  AV 


(2-12) 


where 


BE.m  =  «e  +  max(0,  -ipv^A)g  „^) 


Bs,m  =  +  wiaA:((pVj,zl)^  0) 

BN,m  =  «n  +  rnax{Q,  -[pVyA)^^) 
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^P,m  ^w,m  ^e,m  e,m  ^P^x^\ 


[pVyA)^^  SjjPjn^V 


a 


e,m 


^  r],e,m 

XE-Xp 


a 


w,m 


^  r],w,m 
Xp—X]^ 


a 


n,m 


^  r],n,m 

yN-yp 


a 


s,m 


^  T],s,m 

yp-ys 


(1  ^2}Vw,m.  ”1"  ^3VwW,m}^^^^(^(.P^xV^^w,m> 

{-kiriE,m  +  (1  -  k2)Vp,m  “  max{ipv^r]A) e,m>  o)  + 

(  ^4VE,m.  "I"  (1  ^5)Vp,m.  ^P^xP^w,m> 

(j^4^EE,m  ^5)VE,m  ”1"  ^P^xP‘^')e,m>  ”1” 

(kjTjp^rn  -  (1  -  ks)ris,m  +  k^rjssm)^^^  (^(pVypA) o)  + 
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i-k7riN,m  +  (1  -  ks)rip,m  -  k^rjs^m)  max  (^(pVyr]A)^^,  o)  + 

{~kio'nN,m  +  (1  “  kii)r]p  jn  — 

ki2ris,m)max  (^-(pVypA)^^,  o)  +  -  (1  - 

kii)r]N,m  +  ki2Vp,m)max  (^-(pVypA)^^,  o)  + 


(2-13) 


In  general,  the  discretized  equation  can  be  stated  in  the  following: 


^ P,mVp,m  ^i=l  ^ ^ 


(2-14) 

In  Equation  (3-14),  n  is  the  number  of  faces  surrounding  the  calculating 
volume.  Hi  is  the  value  of  q  passing  through  face  i,  Sr,  is  the  source  term. 

2.3  Multiphase  model 

Since  there  are  two  phases  involved  in  the  system,  the  multiphase  modeling 
is  required.  In  this  part,  two  different  modeling  methods  will  be  introduced. 
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2.3.1  Mixture  Model  Theory 


Continuity  equation: 

^  +  V  ■  (PaveVave)  =  «  (2-15) 

In  Equation  (3-15),  Pave  is  the  average  density  of  the  mixture,  and  Vave  is 
the  average  mass  velocity.  The  average  density  and  mass  velocity  of  the  mixture 
can  be  expressed  as  follows  [16]: 


Pave  Si=l  ^iPi 


(2-16) 


^ave 


_  Id=iaiPiVi 


Pave 


(2-17) 


Where  tti  is  the  volume  fraction  of  phase  i. 


Momentum  equation: 


(Pave^ave)  ”i"  ^  '  ^P ave^ ave^ ave) 


^  '  [Pave  {^^ave  ”i"  ^^ave  ”i"  Pave% 
~^ave  ^  ^  '  ^iPi^r,i  ^r,i 


i=l 
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(2-18) 


Where  F  is  the  body  force.  Have  the  average  viscosity  of  the  mixture 
and  Vy  i  is  relative  velocity  between  phase  I  and  the  average  mass  velocity. 


We  can  define: 


Vsp  =  Vs  -Vp 

where  is  the  velocity  of  the  secondary  phase,  and  Vp  is  the  velocity  of 
the  primary  phase.  There  is  another  parameter,  mass  fraction  of  any  phase  i,  that 
needs  to  be  defined: 

_  ^ 

^  i 

Pave 

Combined  these  terms  results  in  the  following  equation: 


Vr,i  =  V 


sp  Si=l  ^iVpi 


(2-19) 
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2.3.2  Volume  of  Fluid  Model  Theory 


Basically,  the  Volume  of  Fluid  Model  Theory  [16]  solves  the  multiphase  fluid 
dynamics  by  tracking  the  ration  of  each  fluid  in  the  control  volume.  If  Ok  is  the 
volume  fraction  of  kth  phase,  following  expression  can  be  obtained: 

=  1:  The  cell  is  full  of  the  kth  phase 

=  0:  The  cell  is  empty  of  the  kth  phase 

0  <  «/£  <  1:  The  cell  contains  an  interface  between  the  kth  phase  and  other 
phases 

The  continuity  equation  of  the  kth  phase  can  be  stated  in  the  following  form: 

^  (tt/cP/c  +  V  ■  {akPjVk))  =  rhik  -  rhki  (2-20) 

where  rhu^  is  the  mass  transfer  from  phase  I  to  phase  k. 

Since  there  is  only  one  momentum  equation  required  in  this  method,  the 
properties  in  the  transport  equations  are  required  to  be  defined.  Generally,  the 
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properties  can  be  calculated  in  the  following  form: 


5  =  (2-21) 

where  6  is  the  property. 

The  momentum  equation  in  the  Volume  of  Fluid  Model  theory  is  given  in 
equation  2-22: 

^(“aPa^a  +  +  V  '  (ttaP^^a^a  +  = 

-Vp  +  V  ■  (VVa  +  Vv'^’a)  +  {Vv^  +  Vv’^’j^)]  +  {UaP^  + 

«wPjg  +  F 

(2-22) 

In  this  research,  momentum  equation  for  the  algae  [17]  phase  is  added: 

PaVa^  =  PaVag+VaVp  +  ¥ 

(2-23) 


Where  [18][19][20] 
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It  follows  that,  more  resistance  is  added  when  flow  through  porous  media  is 
considered.  In  this  study,  a  resistance  term  is  added  to  the  momentum  equation 
for  the  fluid  dynamics  in  porous  media.  The  resistance  term  can  be  dealt  with  as  a 
negative  effect  source  term  in  the  modeling,  and  it  can  be  expressed  in  the 
following  form  [16]: 

Sj  =  -  (SLi  (^j,iP  \v\Vi) 

(2-25) 
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where  Sj  is  the  source  term  for  the  jth  momentum  equation,  and  |v|is  the 
magnitude  of  the  velocity.  The  first  term  on  the  right  hand  side  is  the  viscous  loss 
term,  and  the  second  term  on  the  right  hand  side  is  the  initial  loss  term.  If  the 
porous  media  is  homogenous,  the  following  source  term  results: 


where  a  is  permeability,  and  C'  is  the  initial  resistance  factor.  Hence,  the 
matrices  D  and  C  can  be  stated  in  the  following  form: 


2.5  Modeling  Surface  Tension 


Although  the  effect  of  surface  tension  effect  has  been  shown  empirically  by 
AVS  to  have  the  same  effect  as  a  4  cm  pressure  head  in  the  capillary  zone, 
mathematical  modeling  is  still  important  to  be  discussed.  Brackbill  [21][22]  and 
others  have  proposed  that  surface  tension  can  be  expressed  as  a  volume  force 
term  that  is,  added  it  to  the  momentum  equation  as  source  term.  The  source  term 
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can  be  stated  in  the  following  form: 


^vol 


aiPiKjVaj+ajPjKjVai 

'j  0.5(pi+pj) 


(2-26) 


Where  CTj  j  is  the  surface  tension  coefficient,  and  Kj  is  the  curvature  of  the 
ith  phase,  which  can  be  defined  in  the  following: 


K 


=  V 


Va, 


(2-27) 


2.6  Pressure  Correction 


Generally,  the  momentum  equation  can  be  discretized  in  the  following  form 
(x  direction): 


B 


P,m+i^P,x,m+l 


II 

1 


i=l 


^  i,m+l^i,x,m+l  +  ^m+1 


fSPm+l\ 

[  Sx  ) 


Since  the  non-linear  equations  of  the  system  are  linearized  in  each  time  step, 
the  entire  iterative  process  occurs  in  two  different  steps.  In  the  outer  iterations, 
the  coefficient  and  source  terms  are  updated,  while  the  linear  systems  are  solved 
with  fixed  coefficients  in  the  inner  iterations.  For  each  outer  iteration,  the 
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equations  can  be  solved  in  the  following  form: 


(2-28) 


Where  V  is  a  temporary  value  of  Due  to  the  reason  that  the 

pressure  term  applied  in  equation  (2-28)  is  from  the  previous  iteration,  the 
velocities  calculated  do  not  satisfy  the  continuity  equation.  Therefore,  the  velocity 
is  required  to  be  corrected  to  enforce  the  continuity  equation.  A  velocity  and 
pressure  correction  can  be  defined  as: 

27^+1  =  Vx  +  V'x 
pin+l  _  pin  _|_  pi 


Define 


(2-29) 
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Since  the  continuity  equation  must  be  satisfied,  the  following  equation 


(discrete  Poisson  equation  for  the  pressure)  can  be  obtained: 

=  0 
p 

The  term  Yi=i^i'^'i  can  be  approximated  to  the  discrete 

Poisson  equation  for  the  pressure  correction  can  be  stated  as  follows: 


S 

.dt 

5xi 

,m+l' 


■‘P 


dp 


dt 


S 

jp  Sxi 


S 

P  (Sp\ 

Idpl 

Sxi 

Bp + i]r=i  v  ) . 

p 

.dt. 

p 


(2-30) 

The  velocity  field  is  updated  when  equation  2-30  is  solved.  This  procedure  is 
known  as  SIMPLEC  [23]  algorithm,  and  the  procedure  can  be  states  as  follows: 

1.  Velocity  field  is  available  at  time  step  m 

2.  Calculate  a  temporary  value  of  velocity  at  time  step  m+1 

3.  Solve  the  Possion  equation  for  the  pressure  correction 
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4.  Solve  the  velocity  correction  equation,  and  calculate  momentum 
equation 

5.  If  converge,  go  to  next  time  step.  If  not,  repeat  the  procedure  until  both 
continuity  and  momentum  equations  converge 


2.7  Solver  for  Discretization  equations 


2.7.1  Gauss-Siedel  method 


In  this  study,  the  solver  adopted  here  is  based  on  the  Gauss-Siedel  method 
[24].  An  introduction  to  this  method  is  provided  next: 


Assume  a  square  system  of  linear  equations  given  by: 


Ax  =  b 


Where: 


-ail 

%2  • 

^In 

-Xi' 

A  = 

<^21 

<^22 

^2n 

,  X  = 

^2 

,  b  = 

^2 

.^nl 

^n2 

^nn. 

_^n. 

bn. 
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The  matrix  A  can  be  decomposed  into  the  sum  of  an  upper  triangular,  and 
lower  triangular  matrix. 

L  +  U 

The  system  equation  can  be  rewritten  as: 

Lx  =  b  —  Ux 

A' can  be  solved  by  using  Gauss-Siedel  method  and  the  value  of  A'from  the 
previous  step.  The  iterative  formulation  can  be  stated  as  follows: 

=  L~^(b  -  Ux^)  (2-29) 

Once  the  calculating  matrix  is  formed,  the  solver  can  have  the  solution  for  every 

iteration. 
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Chapter  3  Web  Tension  Control  System  Design 


Web  tension  is  a  key  factor  to  controlling  fluid  dynamics  in  the  capillary 
zone.  In  this  chapter,  we  develop  a  web  tension  control  system  for  the  HDD 
system.  There  are  several  issues  that  need  to  be  considered  in  the  design  of  the 
HDD  system.  To  begin  with,  the  material  chosen  for  the  capillary  belt  should  be 
able  to  draw  as  much  water  from  the  algae  belt  as  possible.  Second,  the  tension 
range  for  the  belt  needs  to  sufficient  to  achieve  the  operating  goals  of  the 
system;  insufficient  tension  can  cause  the  web  sag,  wrinkle  and  rail,  and  if  the 
tension  is  too  high  it  can  cause  deformation  or  breakage.  The  tension  and  the 
transport  velocity  combine  to  define  the  operating  characteristics  of  the 
capillary  belt  and  need  to  be  controlled  simultaneously  in  order  to  achieve  good 
performance. 

3.1  Web  Transport  System 

3.1.1  Wind/Unwind  Transport  System 

There  are  two  dynamic  behaviors  in  the  transport  process:  lateral  and 
longitudinal.  The  longitudinal  behavior  is  due  to  the  variation  of  tension  and 
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transport  velocity,  while  the  lateral  behavior  is  due  to  the  variation  of  lateral 


displacements  in  the  transport  process.  From  Hook's  law,  the  web  deformation 
is  directly  related  to  tension  and  velocity.  In  a  conventional  wind/unwind 
process,  the  rotational  inertia  on  the  wind  and  unwind  rolls  changes  with  radius 
[40],  and  the  radius  can  be  calculated  by  knowledge  of  the  angular  velocity. 
Before  setting  up  the  model,  there  are  some  assumptions:  (1)  the  web  material 
is  an  ideal  elastic  body,  (2)  there  is  no  slip  age  between  the  roller  and  the  web, 
(3)  web  tension  caused  by  the  unwind  roller  is  negligible,  and  (4)  the  torque  due 
to  coulomb  friction  on  the  unwind  roll  is  negligible. 

Figure  3.1  shows  the  roller-to-roller  control  system  in  a  conventional 
wind/unwind  process  where  the  web  material  is  transferred  from  the  unwind 
roller  to  the  rewind  roller.  According  to  the  mass  conservation  law: 

(3-1) 

Temporal  change  of  mass  in  the  control  volume  given  on  the  left  hand  side 
of  the  equation  are  equated  to  the  difference  between  the  input  and  output 
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mass  of  the  control  volume  given  on  the  right  hand  side  of  the  equation. 


Therefore,  the  following  nonlinear  equation  can  be  obtained: 

dt  \l  +  £jk(t))  l  +  £jk(t)  l  + 


(3-2) 


where 


d  (  i-ijkO-^  _  1 

dt  \l+£jk(t))  l+£jk(t)  dt  {i+£jj^(t)f  dt 

And  p  Density  of  the  web 

V  Reference  velocity  of  the  web 

Lji^  The  length  of  the  web  between  the  rollers,  where  j  represents  the 
unwind  roller  and  k  represents  the  rewind  roller. 

£ji^  The  strain  of  the  web  between  rollers 

Substitute  equation  (3-2)  into  equation  (3-1),  we  have  the  following  equation: 

(l  +  EjM)  ^  -  LjM  =  (l  +  EjM)  -  vM) 


41 


(3-3) 


According  to  Hook's  Law,  the  relationship  between  strain  of  the  web  and  the 
tension  can  be  described  as  follows: 


T  =  AEs 


Therefore,  following  equation  can  be  obtained: 

V  AE  /  dt  AE  dt  V  AE  /  \  J 

(3-4) 

(3-5) 


Here,  the  spring  constant  =  —  ,  and  the  above  equation  can  be  rewritten  as 

^;7c 


follows: 


(3-6) 


In  the  wind  process,  the  radius  of  the  rewind  roll  changes  continuously,  and  the 
rotational  inertia  varies  with  the  radius.  The  dynamics  of  the  unwind  radius  is: 
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(3-7) 


hvj(t') 

2nRj(t) 


And  the  dynamics  of  the  rewind  radius  is: 


RM 


hVkit) 

2nRk{t) 


(3-8) 


Where 

R  Radius  of  roller 

h  The  thickness  of  the  web 

V  The  velocity 


The  rotational  inertia  of  the  rolls  can  be  controlled  by  the  motor  torque  [40]  [42], 

and  the  state  equation  can  be  described  as  follows: 

d  _  dtfi(t) 

dt  Rj(t)  dt 

Vj(t)  Vj(t) 

=  -^J  +  TjMRjit)  =  KjUj  -Dj-^  +  TjMRjit) 


(3-9) 


Where 
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Tj(N  ■  m)  Output  torque  of  unwind  roll 


Kj(N  ■  Unwind  motor  torque  constant 

S6C 

Dj(N  ■  m - )  Friction  coefficient  of  unwind  roll 

^  ^  rad-' 


Since  the  cylinder  is  composed  of  the  web  and  metal,  the  rotational  inertia  of 


the  cylinder  can  be  described  as  follow: 


J j  J jm 


phy^Tl(R:{t) 


(3-10) 


Where  is  the  inertia  of  the  metal  roll,  is  the  radius  of  the  metal  roll. 


and  is  the  width  of  the  web.  Therefore,  the  differential  equation  of  the 


rotational  inertia  will  be: 


d]j(t)  _ 
dt 


4hK^Rj(t) 

2tt 


(3-11) 


Where  and  the  following  equation  can  be  obtained: 


u(t)  = 


^hKyyRf  (t)  vf  (t)  )  hvf  (t) 


Rj{t)KjUj  -  DjVjit)  + 


+ 


2n  J  2nRj'  (t) 


(3-12) 


Similarly,  the  differential  equation  for  the  rotational  inertia  of  the  rewind  roll 
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can  be  described  as  follow: 


4(t) 


1 

Jkm  +  -  Rtm] 


4hK^Rlit)vl(t)]  ^  hvlit) 
2n  j  2nRl(t) 


(3-13) 


Fig.  3.1  Web  tension  control  system 
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Fig.3.2  Simplified  web  tension  system 
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Mi 


Ri  :  Tension  Regulated  Motor 
Rj  :  Tension  Regulated  Slave  Motor 
Rk :  Speed  Regulated  Head  Motor 


Head 

Pulley 

(Speed 

Pacer) 


3.1.2  HDD  web  transport  system 


Figure  3.4  shows  the  HDD  web  transport  system,  including  four  rollers  and 
tension  regulators  Ri  and  Rj  that  are  used  to  control  web  tension;  Rk  is  a  velocity 
regulator,  which  is  used  to  control  transport  velocity;  Rm  is  the  displacement 
regulator,  which  is  used  to  maintain  the  tension  between  tail  and  head  pulleys. 
In  this  system,  the  radius  of  rollers  can  be  considered  as  constant  since  there  is 
no  wind/unwind  process  involved.  The  governing  equations  can  be  adjusted  as 
follows: 
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.  AE  J  dt 


AE 


dt 


(3-14) 

^  -  RjOJjCt))  +  -  RjOJjCt)) 

(3-15) 


Here,  the  spring  constant  =  —  ,  and  the  above  equation  can  be  rewritten  as 

^;7c 


follows: 


^  =  %  {R,co,(t}  -  Rjcojit))  +  {R,co,(t}  -  RjCOj(t)) 

(3-16) 


Using  the  same  method. 


^  =  Ki,  (Rjcojit)  -  Ricotit})  +  '!J^[Rjco,(t}  -  R,co,(t}) 


(3-17) 


=  RuiiRif^iit)  -  Ru<iu(t))  -  R^oiuit))  + 

_l_  ^fct(t)\  dLjf^i(t') 

V  ifeift)/  dt 


(3-18) 
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Since 


^ki 


=  Lit  +  V (X;  -  Xfc)2  +  (y,(t)  -  yfc)2  -  (fi,  -  +  7 (X;  -  +  (y,(t)  -  yt)2  -  (R,  - 


(xi  -  +  (yi(t)  -  y^Y  -  (Ri  -  R^Y  +  (x,  -  +  (yijt)  -  -  (R;  -  -  c] 

27 -  Xfc)2  +  (y,(t)  -  yfc)2  -  (R,  -  Rfc)y  (x,  -  +  (y,(t)  -  yt)^  -  (R,  -  1 


=> 


rfLfci(t)  yi(t)-yfc  .  ,  yi(t)-yt 

dt  7“  +  (.yi(t)-yk)  7^^  +  (yiW->^ 


a  +  (yi(0  -  yfc)^  +  fc  +  (yi(0  -  yt)^  -  cy]  ^ 

.  27“  +  (y;  (0  -  yfc) +  (y;  (0  -  ytY  /  j 


1 

2(a  +  (y,(t)  -  yfc)2)(&  +  (y,(t)  -  yt)^) 


|(4y,(t)  -  2yfc  -  2yt)^ja  +  {yi{t)-y^y^b  +  {yi{t)  -  ytYyiit) 

-  (a  +  (y,(t)  -  yfc)2  +  fo  +  (y,(t)  -  yt)^ 

((yi(0  -  yfc)7^  +  CyiCO  ~  yt)^ 

(  7“  +  (yi(0-yfc)^ 

,  (yi(t)  -  yt)7«  +  (yi(0  -  yfc)^] . 

+ - r  = - f  yi  (t)  f 


(3-19) 


Where 
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a  =  (xj  -  Xfc)2  -  (Ri  -  R^y 


XtY  -  (Ri  -  Rty 
c  =  (xt  -  x^y  +  (yt  -  yjc)^ 
yi(t)  = 

The  rotational  inertia  of  the  rolls  can  be  controlled  by  the  motor  torque,  and  the 
state  equation  can  be  described  as  follows: 
d 

—Jj(x)j(t)  =  -Tj  -  DjOijit)  +  (Tjkit)  -  Tij{t))Rj 

=  KjUj  -  Djoijit)  +  -  TtjityRj 


(3-21) 
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d 


-Tj  -  DiOiiit)  +  (Jijit)  - 


=  KiUi-  DiOiiit)  +  (Tijit)  -  T^iimRi 


Where 


T(iV  ■  m)  Output  torque  of  roll 


K{N  •  Motor  torque  constant 


(3-22) 


(3-23) 


D(N  -m - )  Friction  coefficient  of  roll 

rad'' 


3.2  Control  Theory  and  Controller  Design 


From  the  mathematical  model  obtained  in  the  previous  section,  we  can  see 


that  the  system  is  highly  nonlinear.  In  this  study,  we  develop  a  controller  based  on 
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sliding  mode  control  (SMC),  which  is  a  well  known  control  method  providing  a 


robust  controller  for  nonlinear  dynamic  systems  [43],  However,  chattering 
phenomena  caused  by  discontinuous  control  will  take  place  while  the  system  is 
operated  near  the  sliding  surface.  One  of  the  most  common  solutions  to  this 
problem  is  to  introduce  a  boundary  layer  on  the  sliding  surface  [44],  The  details 
of  the  design  of  the  controller  are  introduced  in  the  following  paragraphs. 

3.2.1  Sliding  Mode  Control 

Sliding  mode  control  is  a  form  of  variable  structure  control  (VSC).  It  is  a 
nonlinear  control  method  that  changes  the  dynamics  of  a  nonlinear  system  by 
applying  a  high-frequency  switching  control  function.  There  are  two  modes;  the 
approaching  mode  and  the  sliding  mode.  In  the  approaching  mode,  the 
trajectory  is  directed  toward  the  sliding  surface  in  limited  time  ,  while  in  the 
sliding  mode  the  trajectory  is  confined  to  the  sliding  surface  where  is  moves  and 
toward  the  equilibrium  point.  The  surface  function  (S)  can  be  described  as 
follows: 

S  =  e(t)  +  ae(t)  +  /?  J  e  dt 
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(3-24) 


where  e(t)  G  is  the  tracking  error,  a  and  |3  are  constants.  Let's  assume  an 
n-th  order  nonlinear  system  that  is  linear  in  the  control  input  (u) 

x(t)  =  /(x,  t)  -I-  5(x,  t)u 

where  x  —  [x, x',  ...,x”“^],  and  =  [x^.x'^, ... ,x^~^]  is  the  desired  value. 
Hence,  the  tracking  error  is: 

e  =  X  —  Xd 

The  sliding  surface  can  be  described  as: 

S{e,t)  =  CQje  (t)dt  +  e’'(t)  +  C„e”(t)  (3-25) 

where  Cq,  C^,  ... ,  are  positive,  p  and  q  are  positive  odd  integers,  and  n  is  the 
relative  degree,  r  is  positive  integer. 

Since  the  sliding  surface  is  the  index  of  tracking  error,  it  can  be  the  index  of 
tracking  performance.  Therefore,  the  control  logic  can  be  described  as  follow: 

<  0,when  S(e,t)  >  0 
|5(e,  t)  >  0,whenS(e,t)  <  0 
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=>  55  <  0 


Since  there  exists  chattering  in  SMC,  the  sign  function  is  replaced  by 
saturation  function  in  this  study  in  order  to  ease  the  chattering  phenomena. 
Therefore,  5  can  be  described  as  5  =  —fj.sat(S),  where  fj.  \s  a  control  gain,  and 
sat(5)  is  the  saturation  function. 

(  1  S>d 

sat(S,  d)  -  |/(5,  d)  -d  <  S  <  d 

(  -1  S<-d 

Then  using  the  control  logic  as  define,  we  can  define  a  candidate  Lyapunov 
function  as: 

Ly(5)  =  “‘5’^  >0,  and  Ly(S)  =  55  <  0 

Therefore,  Ly  is  uniformly  stable,  S  is  bounded,  and  5  G  L2.  Since,  each 
term  of  S  is  bounded,  and  S  E  L2,  5  is  bounded.  From  Barbalat's  lemma,  the 
following  conclusion  can  be  reached: 

limt^ool^l  =  lim|x-Xdl  =0 

t->CO 


54 


3.2.2  Sliding  Mode  Controller 


In  this  study,  the  tracking  error  is  defined  as  follow: 

=  ^;7c(t)  “  Td 

The  differential  of  the  tracking  error  of  the  web  tension  is: 

%(t)  =  ^;7c(0  =  ^  +  - 

/?yCOy(t)) 

(3-26) 

%(0  =  Tj„it)  = 

^fert)fc(t)T'yfc(t)  +  i?fcajfc(t)'ryfc(t)  — i?y  (i)y  (t)T'yfc(t)— i?ya)y  (t)'ryfc(t)  TT"  ^  f)  • 

- : - r  - 

^jk 

Rjthjit}) 

(3-27) 

Now,  we  can  define  the  sliding  function. 

Srj.it)  =  er.^t)  +  +  Prj,  I  erj,  dt 

STj.dt)  =  erj.it)  +  ar.^er.^t)  +  pT,^er.^{t)  =  -iir.^sat{Sr.^) 
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->  - 

Rkd)kit}Tjkit)  +  Rka)k(t)tjkit}  -  Rjd}j(t)Tjkit)  -  Rja)j{t)fjk{t) 


+  K(Rkd}kit)  -  Rjd)jit))  +  aj  fjk(t) 


(  ^jTjk(^)\  .  .  . 


KR,  +  ^  -  rj 


->  - 


{Rk^k  -  Dk(^kii)  +  (T’feiCO  -  7’^fe(t))i?fe}  - 


{  ^  hk 


[KjUj  -  Dj(x)j{t)  +  O' jkA)  -  7’y(t))i?J  +  (a^.^  + 


Li. 


^O  jkiO 


)  (^k^^kCO  - 


-  Td]  =  -tir.O^tOTjA 


Hence,  uj  can  be  chosen  as  follows: 
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Uj  = 


-  (TjM  -  Tijit))  Rj  + 


Ji 


KR],+ 


^]k 


RjTjUt) 

KRi+ 

^  ^jk 


Jk 


\j^k^k 


D,a),(t)  +  -  r;,(t))i?  J  +  ( arj,  +  [k 


^jk 


RjCDjCt)^  +  ^^(^Rk(i)k(.i)  —  Rj(i)j(.t)^  +  PTj^[Tjk(.t) -Ta\  + 


fiTj^satdSrJ 


jk^ 


(3-28) 


Similarly,  Ui  can  be  chosen  as  follows: 
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+ 


Ji 


KRk-^  7 


kj  ) 


KRi  + 


RJjjit) 


J, 


\KjUj  —  DjCxijit) 


+  -  Tij(t))Rj\ 

R,a)i{t)  RiTuity 


l.  I  I 


K  [Rja)j(.t)  -  Ria)i(.t)') 
+  PTjTijit)-Ta] 


+  ^^T^Jsat(ST^J) 


(3-29) 


Define  the  velocity  error  function 


Rk^kO'^  Rk^kd  ^  Rk^kO'^ 


The  sliding  function  can  be  defined  as  follows: 


Sv(,t)  =  ey{t)  +  pyj  Cy  (t)dt 

^  ‘^vCO  ^vCO  ”h  Py6y  (t)  fiySCLt(^Sy^ 
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=>  5v(t)  =  -  DfcCOfc(t)  +  (rfcj(t)  -  7;fc(t))i?fc]  +  - 


^k^kd)  liySCLt(Sy^ 


(3-30) 


Hence,  u^.  can  be  chosen  as: 

U/c  = 

-  (r^jCO  -  r^fe(t))  ^  {py{Rk(^k(i)  -  Rk^^kd)  + 


(3-31) 

Since  is  adopted  to  control  the  displacement  of  displacement  regulator, 
and  the  only  concern  is  that  the  tension  between  ith  and  kth  roller  has  to  be 
bigger  than  zero,  the  controller  can  be  chosen  as  follows: 

(  0.01  T„i  <  O.SN 

Um=^,  0  ,  O.SN<ni<lN 

™  (-0.01  T„i  >  IN 
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(3-32) 


Chapter  4  Experimental  Methods 


One  of  the  main  objectives  of  the  study  is  to  establish  analytical  model  for  the 
hydraulic  system  of  the  HDD  process  in  both  the  gravity  well  and  capillary  zone. 
This  requires  some  characteristics  studied.  First  of  all,  it  is  important  to  know  if  the 
fluid  is  laminar  or  turbulent.  If  the  fluid  dynamics  is  laminar,  then  the  general  form 
of  Darcy's  Law  can  be  adopted  in  the  modeling.  If  it  is  not,  a  turbulent  flow  model 
is  required.  Second,  the  permeability  in  the  gravity  well  will  change  with  the 
volume  of  fluid,  so  it  is  necessary  to  quantify  this  characteristic  in  the  model.  Third, 
temperature  on  the  fluid  dynamics  is  a  concern.  Since  the  algae  are  alive  during 
the  dewatering  process,  the  temperature  may  change  the  fluid  dynamics 
properties  of  the  system.  Hence,  the  temperature  effect  needs  further  study.  In 
this  study,  there  are  a  series  of  experiments  that  are  done  for  testing  these 
characteristics.  The  experimental  method  and  results  will  be  provided  in  this 
chapter. 

4.1  Permeability  test 

From  previous  chapter,  the  governing  equations  can  be  formulated  in  any 
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arbitrary  geometry.  However,  some  key  factors  are  still  unknown  in  the  model.  The 


main  purpose  of  this  chapter  is  to  define  and  quantify  them  to  complete  the  CFD 
modeling.  Figure  4.1  provides  the  basic  idea  of  the  experiments,  with  the  objective 
of  determining  the  correct  fluid  dynamic  model  to  use.  In  particular,  the 
experiments  are  designed  based  on  Darcy's  Law.  As  we  know,  Darcy's  Law  was 
formulated  based  on  the  result  of  experiments  and  derived  from  the  Navier-Stokes 
equations  via  homogenization,  and  it  has  been  used  to  describe  the  fluid  dynamic 
of  flow  through  porous  media  in  a  wide  array  of  applications.  In  this  experiment, 
Darcy's  Law  can  be  simplified  in  the  following  equation  if  the  driving  force  and  the 
temperature  are  constants. 


-kA  AP 
II  L 


(4-1) 


Based  on  this  idea,  the  permeability  of  the  filters  and  algae  can  be  estimated 
by  simply  measuring  the  flow  rate.  In  figure  4.1,  there  are  several  control  modules 
that  regulate  the  input  operating  conditions.  The  basic  idea  of  the  experiment  is  to 
feed  algae  solution  into  a  testing  tube,  and  deposit  the  algae  on  the  filter  to  reach 
the  desired  amount  of  wet  algae.  Then,  water  is  fed  into  the  filter/algae  at  a 
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desired  pressure  head,  where  the  pressure  head  control  module  controls  the 


solution  temperature  to  be  constant.  In  these  series  of  experiments,  temperature. 


pressure  head,  and  thickness  of  algae  are  our  main  concerns. 


Fig.  4.1  Experiment  method 


A  filter  with  15  urn  pore  size  and  9%  porosity  was  chosen  for  the  testing.  First  of  all. 


experiments  to  determine  the  permeability  of  15-9  filter  were  conducted.  Figure 


4.2  shows  the  measurement  of  flow  rate  during  the  testing  process.  As  you  can  see. 


the  flow  rate  is  very  constant,  and  the  Reynolds  number  in  this  case  is  around 
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2252.  Therefore,  the  flow  is  laminar,  and  Darcy's  Law  can  be  applied  in  the  tests. 


Figure  4.3  demonstrates  the  measurement  of  flow  rate  for  different  pressure 
heads.  Figure  4.4  shows  the  correlation  between  resistance  (pressure  head  divided 
by  flow  rate)  and  pressure  head,  and  it  shows  that  resistance  is  not  changing.  The 
temperature  was  controlled  at  75°F,  and  the  total  amount  of  algae  deposited  on 
the  filter  is  around  0.105  g.  Again,  the  flow  rate  remains  constant  for  each  pressure 
head,  and  it  is  highly  linear  to  the  pressure  head.  Hence,  the  experiments  are 
consistent  with  Darcy's  Law,  pressure  drop  positively  linearly  correlated  to  flow 

rate. 


Since  flow  rate  is  related  to  thickness  of  the  porous  media,  the  measurement 
of  algae  thickness  is  required.  Because  the  thickness  of  wet  algae  changes  with 
moisture  content,  the  effect  of  changes  in  moisture  content  on  algae  thickness 
must  be  estimated.  In  this  series  of  experiments,  both  the  flow  in  the  gravity  well 
and  the  capillary  belt  are  investigated.  To  begin  with,  the  gravity  well  environment 
was  simulated.  The  idea  is  to  feed  the  algae  solution  into  the  simulated 
environment,  and  deposit  algae  on  the  filter.  Then,  water  is  fed  into  simulated 
gravity  well,  and  the  thickness  of  the  wet  algae  is  observed.  As  we  know,  the  algae 
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solution  keeps  feeding  into  the  gravity  well  in  AVS  system,  and  given  the  condition 


that  resistance  grows  rapidly,  the  moisture  content  of  the  algae  deposited  on  the 
belt  will  not  change  significantly.  From  experimental  observations,  the  thickness  of 
the  algae  didn't  change  with  time  passing  or  with  changes  in  pressure  head  (fig. 
4.5).  Hence  there  is  no  need  to  explore  the  thickness  module  in  gravity  estimation 
model  any  further.  In  the  second  part,  the  phenomenon  in  the  capillary  zone  is 
explored.  Since  the  algae  belt  in  capillary  zone  is  oriented  with  an  upward  slope, 
the  algae  needs  to  reach  a  certain  moisture  content  to  stay  on  the  belt.  According 
to  the  experimental  data,  the  moisture  content  has  to  be  less  than  99%,  and  it  has 
to  be  in  the  form  of  sludge.  The  procedure  is  to  use  a  vacuum  filtration  system  to 
make  samples,  and  use  laser  profilometer  to  measure  the  thickness  of  the  samples. 
Usually,  vacuum  filtration  is  faster  than  a  simple  gravity  filtration  because  the 
vacuum  system  creates  a  large  driving  force  between  the  media.  The  vacuum 
filtration  system  setup  is  shown  in  figure  4.6.  A  laser  profilometer  (Fig.  4.7)  is  a 
non-contact  measurement  device  that  can  measure  surface  roughness  and  film 
thickness  if  a  ground  reference  is  given.  The  calculation  of  roughness  is  given  in  the 
following  equation: 
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R 


rms 


(4-2) 


Where  Rrms  is  the  surface  roughness  on  the  basis  of  the  root  mean  square 
deviations,  where  yi  is  the  difference  in  height  between  i-th  point  and  mean  value. 
The  surface  roughness  of  the  algae  sample  is  around  2.79  urn  while  the  thickness 
of  wet  algae  is  around  327  urn.  Hence,  the  assumption  that  the  thickness  of  the 
algae  sludge  is  very  uniform  is  valid.  Figure  4.8  is  the  measurement  of  wet  algae 
thickness  for  different  moisture  contents.  As  the  figure  shows,  the  thickness  of  wet 
algae  in  our  range  of  interesting  doesn't  change  too  much  although  it  shrinks  very 
fast  thereafter.  From  these  experiments  we  conclude  that  it  is  reasonable  to 
assume  that  thickness  of  the  wet  algae  is  approximately  constant  in  both  the 
gravity  and  capillary  zone.  In  order  to  properly  include  thickness  into  the  model, 
the  resistance  part  of  the  model  needs  to  be  replaced  by  another  parameter, 
which  should  be  a  function  of  the  depositing  area,  viscosity,  concentration  of  the 
algae  solution  and  permeability.  Before  developing  this  model,  there  are  some 
parameters  that  require  further  investigation.  First  of  all,  temperature  is  known  as 
a  key  factor  that  can  influence  flow  rate  through  changes  in  the  viscosity  of  water 
and  perhaps  characteristics  of  the  algae.  The  experimental  setup  is  just  like  Fig.  4.1. 
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Figures  4.9  and  4.10  show  the  results  of  the  temperature  tests.  The  flow  rate  was 


monitored  at  75,  88,  96,  and  106°F,  the  pressure  head  is  controlled  at  13.25  inches 
of  water,  and  the  total  amount  of  algae  was  controlled  at  the  same  level  for  each 
test.  The  results  show  that  the  flow  rate  doesn't  change  too  much  across  the 
temperature  range.  However,  we  know  that  the  dynamic  viscosity  of  water 
changes  with  temperature.  The  dynamic  viscosity  at  75°F  is  around  8.9091E-4 
kg/m»s,  while  it  is  6.2232E-4  kg/m»s,  which  is  more  than  a  20%  difference,  which 
means  the  permeability  of  algae  should  be  changing  with  temperature.  Either  the 
change  in  resistance  is  too  small  to  be  determined  by  our  experiments  or  the  algae 
are  somehow  be  able  to  adjust  to  maintain  a  constant  resistance.  Permeability 
tests  for  different  algae  amounts  were  done  next.  There  were  19  samples  made  for 
these  series  of  experiments.  The  weights  of  the  algae  were  from  0.007  to  0.075  g, 
and  were  uniformly  deposited  on  a  filter,  11.39  cm^  pass  through  area.  Figure  4.11 
shows  the  relationship  between  the  total  amount  of  algae  and  resistance.  It  looks 
as  if  they  are  linearly  correlated  across  the  testing  range.  As  we  know,  resistance 
should  be  a  function  of  viscosity,  permeability,  and  thickness  of  the  porous  media. 
In  this  series  of  experiments,  since  the  temperature  was  maintained  at  75°F,  and. 
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theoretically,  permeability  of  porous  media  is  constant,  the  resistance  is  only  the 


function  of  thickness  of  porous  media.  From  the  results  shown  in  fig.  4.11,  the 


resistance  is  linearly  related  to  the  total  amount  of  algae  deposited  on  pass 


through  area.  This  is  consistent  with  Darcy's  Law,  and  we  conclude  that  the 


permeability  of  the  wet  algae  is  constant  throughout  the  testing  range. 
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Fig.  4.2  Permeability  test  for  15-9  filter 
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Fig.  4.3  Flow  rate  v.s  Pressure  head  for  15-9  filter 
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Fig.  4.4  Resistance  v.s.  Pressure  head  for  15-9  filter 
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Fig.  4.5  Observation  of  thickness  changing  over  pressure  head 


Fig.  4.6  Vacuum  filtration  system  setup 
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Fig.  4.7  Laser  Profilometer 


Fig.  4.8  Thickness  of  wet  algae  v.s.  moisture  content 
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Fig.  4.9  Temperature  test  for  15-9  filter 
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Fig.  4.10  Resistance  v.s.  Temperature  for  15-9  filter 
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Fig.  4.11  Resistance  v.s.  Algae  weight  over  deposit  area 

4.2  Free  Flow  Tests  and  CFD  Mdeling 

In  order  to  test  the  CFD  model,  free  flow  tests  were  run.  The  idea  of  the 
experiments  is  to  feed  solution  into  the  testing  device,  and  monitor  the  flow  rate 
and  free  surface  during  the  entire  processes.  The  experimental  data  will  be  used  to 
compare  the  simulation  results  based  on  the  CFD  model  discussed  in  the  previous 
chapter.  The  testing  device  in  this  experiment  will  be  simple  and  easy  to  establish 
the  calculation  grids  for  the  CFD  model.  Therefore,  the  simulation  will  be  less 
complex  than  the  real  HDD  system.  Before  testing  the  real  fluid  behavior  in  the 
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simple  HDD  system  model,  water  will  be  used  for  several  different  filters.  Table  4.1 


is  the  permeability  table  obtaining  by  the  permeability  testing.  If  the  fluid  used  in 
the  tests  is  water,  the  following  analytical  equation  can  be  obtained  according  to 
Darcy's  Law. 


A 


dh 

dt 


-kA  AP 
[i  L 


-kA  pg 

==>  h  =  /ioexp( - — ) 

U,  L/ 


(4-3) 


Table  4.1  Permeability  divided  by  L  of  each  filter 


Filter 

Permeability/L 

1.50912E-08 

6_5 

5.85138E-08 

7_2 

1.34753E-08 

18_13 

1.04091E-07 

25_20 

2.84981E-07 

6_3 

6.98285E-09 

In  this  case,  the  control  volume  is  not  constant  because  the  flow  surface 
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keeps  changing  during  the  test.  Hence,  a  dynamic  control  volume  has  to  be 


adopted.  The  CFD  modeling  method  is  demonstrated: 


Since  the  boundary  is  changing  in  this  case,  and  the  only  changing  boundary 
is  the  free  water  surface,  the  spring-based  smoothing  method  can  be  used.  The 
edges  between  any  two  calculating  nodes  can  be  modeled  as  a  spring  using  the 
following: 

Fi  =  S"  kij  (AXj  -  AXi)  (4-4) 

Where  Axj  and  Axj  are  the  displacements  of  node  i  and  neighboring  nodes 
j,  kij  is  the  spring  constant  between  nodes  i  and  j.  The  spring  constant  can  be 
approximated  using  the  following  form: 


(4-5) 


According  to  Hooke's  law,  the  force  generated  between  any  two  nodes  is 
proportional  to  the  spring  constant,  and  the  net  force  will  be  zero  when  it  reaches 
equilibrium.  Hence,  the  iterative  equation  can  be  represented  in  the  following 
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form: 


n+l 


l^kijAxi^ 


(4-6) 


After  updating  the  boundary  nodes,  displacements  can  be  calculated,  and  all 
the  interior  nodes  can  be  obtained  by  solving  Equation  4-6.  Hence,  the  position  at 
the  next  time  step  can  be  determined  from  the  following  equation: 


->  m+l 


=  X; 


m 


+  AXj 


m 


(4-7) 


The  dynamic  meshes  are  shown  in  the  figure  4.12,  and  the  calculation  and 
the  modeling  method  are  described  in  chapter  2.  The  simulation  results  are  shown 
in  figs.  4.13  to  4.18. 


Comparing  the  free  water  flow  experiments  and  the  simulation  results  from 
figs.  4.13  to  fig.  4.18,  we  observe  that  the  three  curves  almost  match  perfectly  in 
each  case.  The  results  show  that  the  CFD  model  can  predict  the  fluid  dynamics  in 
free  water  flow  tests  very  well.  Now,  we  need  to  consider  the  algae  flow.  As  we 
know,  as  more  fluid  passes  through  filter,  resistance  will  increase.  Therefore,  the 
permeability  in  free  algae  fluid  test  is  dynamics  and  as  a  result,  the  governing 


75 


equations  are  nonlinear,  and  obtaining  an  analytical  solution  is  not  possible. 


However,  the  algae  flow  in  HDD  system  can  be  considered  homogeneous  since  it  is 
well  mixed.  Hence,  the  assumption  that  algae  are  deposited  uniformly  on  the  belt, 
and  the  amount  of  algae  deposited  on  the  belt  is  only  related  to  the  total  flow 
passing  through  the  filter.  In  order  to  substantiate  this  assumption,  several 
experiments  are  conducted.  The  idea  is  to  deposit  algae  on  the  filter,  and  use  the 
laser  profilometer  to  measure  the  roughness  and  thickness  of  the  deposited  algae. 
A  laser  profilometer  is  a  non-contact  measurement  device  that  can  measure  the 
roughness  of  a  surface,  and  the  thickness  of  a  film  if  a  ground  reference  is  given. 
As  previously  discussed,  the  assumption  that  the  thickness  of  the  algae  sludge  is 
uniform  is  valid  since  experimental  data  yields  a  roughness  of  around  2.79  urn  with 
a  thickness  of  wet  algae  of  around  327  urn. 
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Fig.  4.12  Dynamic  meshes 

Moreover,  from  measurements,  we  know  that  the  average  size  of  algae  is 
around  15~18  micrometer.  Therefore,  the  pore  size  of  filter  should  not  be  greater 
than  15  urn  to  prevent  the  algae  from  flowing  through  the  filter.  In  this  test,  there 
are  only  four  filters  that  meet  this  requirement.  Therefore  only  the  1-2  filter  and 
the  7-2  filter  were  tested  to  validate  the  CFD  modeling.  Figs.  4.19~4.26 
demonstrate  the  relationship  between  the  experimental  data  and  CFD  modeling 
results.  As  can  be  seen,  CFD  model  is  good  at  predicting  the  fluid  dynamic  for  both 
filters.  According  to  the  data,  the  free  water  surface  height  can  be  predicted  very 
well  for  the  different  filters.  Therefore,  we  may  conclude  that  CFD  method  used  in 
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this  study  can  be  a  good  predictor  for  HDD  system. 


1-2  filter 


Fig.4.13  Simulated  v.s.  measured  changing  heights  of  free  water  test  (1-2  filter) 


6-5  filter 


Fig.  4.14  Simulated  v.s.  measured  changing  heights  of  free  water  test  (6-5  filter) 
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7-2  filter 


Fig.  4.15  Simulated  v.s.  measured  changing  heights  of  free  water  test  (7-2  filter) 


18-13  filter 


Fig.  4.16  Simulated  v.s.  measured  changing  heights  of  free  water  test  18-13  filter) 
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25-20  filter 


Fig.  4.17  Simulated  v.s.  measured  changing  heights  of  free  water  test  (25-20 

filter) 


6-3  filter 


Fig.  4.18  Simulated  v.s.  measured  changing  heights  of  free  water  test  (6-3  filter) 
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Fig.  4.19  Simulated  results  and  measured  data  (O.Olg/L,  1-2  filter) 


Fig.  4.20  Simulated  results  and  measured  data  (0.02g/L,  1-2  filter) 
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Fig.  4.21  Simulated  results  and  measured  data  (0.04g/L,  1-2  filter) 


Fig.  4.22  Simulated  results  and  measured  data  (0.08g/L,  1-2  filter) 
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Fig.  4.23  Simulated  results  and  measured  data  (O.Olg/L,  7-2  filter) 
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Fig.  4.25  Simulated  results  and  measured  data  (0.04g/L,  7-2  filter) 


Fig.  4.26  Simulated  results  and  measured  data  (0.08g/L,  7-2  filter) 
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Chapter  5  Results  and  Discussion 


As  discussed  earlier  in  this  thesis,  the  main  purpose  of  the  study  is  to  develop 
a  predictive  model  of  the  HDD  process  that  can  be  used  in  an  optimal  control 
system  for  producing  high  quality  dried  algae  for  fuel  manufacturing  facilities.  The 
strategy  is  to  take  advantage  of  CFD  modeling  to  predict  the  fluid  dynamics  in 
different  operating  conditions  and  to  develop  a  lookup  table  that  can  be  used  for 
monitoring  and  control  of  the  HDD  system.  There  are  several  control  modules 
including  web  tension  control,  input  (algae  feedstock)  control,  and  transport 
control  in  the  design.  The  web  tension  control  module  and  CFD  simulation  method 
have  been  demonstrated  in  the  earlier  chapters,  and  the  simulation  results  will  be 
discussed  in  this  part.  Once  the  lookup  table  is  complete,  the  optimal  control 
system  for  the  HDD  system  can  be  deployed.  The  simulation  process  and  the 
results  of  simulation  study  will  also  be  discussed  in  this  chapter. 

5.1  CFD  Simulation  and  Results 

Like  the  discussions  in  Chapters  2  and  4,  CFD  modeling  is  an  important  part  of 
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this  study.  CFD  modeling  and  simulation  can  save  time  and  effort  by  eliminating 
unnecessary  trial-and-error  investigations.  In  our  research,  we  use  CFD  to  model 
the  system  in  different  operating  conditions,  and  build  up  a  lookup  table  that  can 
be  used  for  the  control  of  the  HDD  process.  The  details  of  the  method  and  theory 
are  given  elsewhere  in  this  thesis,  the  simulation  results  are  to  be  given  in  this 
chapter. 

5.1.1  Grid  testing 

The  objective  is  to  develop  spatial  and  temporal  grids  that  can  avoid  excess 
calculations  and  reduce  the  computational  time  and  resources  required.  This 
objective  must  be  balanced  by  the  need  to  manage  the  errors  caused  from  an 
inadequate  computational  grid.  Therefore,  choosing  an  appropriate  number  of 
calculating  grids  is  an  important  aspect  of  the  computational  development.  In 
this  regard,  the  time  step  and  spatial  dimensions  both  need  to  be  considered. 
The  time  grid  is  set  up  from  At  =  1/2  to  At  =  1/64,  and  the  spatial  grids  are  from 
20x20  ~  320x320.  Because  the  flow  is  laminar  and  the  geometry  is  not 
complicated,  the  CFD  simulations  are  all  chosen  to  be  two-dimensional.  Figs. 


86 


Q(nT"/s) 


5.1~5.3  are  the  performance  between  different  calculating  grids  for  different 
concentrations  of  algae  in  the  fluid  testing.  As  can  be  observed,  if  the  spatial 
grid  is  160x160  and  ^^t  =  1/32,  this  is  sufficient  for  our  work  since  there  is  less 
than  a  2%  difference  with  the  320x320  grid  and  2lt  =  1/64.  However,  the 
calculation  time  is  reduced  by  1/8  when  compared  to  the  finer  grid.  Hence,  we 
conclude  that  1/32x160x160  grid  is  sufficient  for  this  study. 


Fig  5.1  Grid  test  for  Ig/L  concentration  fluid  (flow  rate) 
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Fig  5.2  Grid  test  for  5g/L  concentration  fluid  (flow  rate) 
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Fig  5.3  Grid  test  for  20g/L  concentration  fluid  (flow  rate) 


S.1.2  Lookup  table 


The  simulation  grid  is  chosen  as  1/32x160x160.  Tables  5.1~5.4  are  the 
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lookup  tables  for  gravity  draining  flow  rate  (GFR),  capillary  draining  flow  rate 


(CFR),  mass  flow  rate  of  the  dried  algae  output,  and  the  remaining  water  at  the 
exit  of  the  capillary  zone  for  different  concentrations  of  the  input  solution  and 
different  web  velocities.  If  an  ideal  heater  is  designed  as  in  fig.  5.4,  a  heating 
condition  table  can  be  developed.  Tables  5.5~5.7  are  the  tables  for  the  heating 
temperature,  required  inlet  air  flow,  and  required  energy  for  producing  1  lb  of 
appropriately  dried  algae.  As  these  tables  show,  it  seems  that  the  concentration 
of  the  input  solution  and  the  web  velocity  are  two  key  parameters  for  the  yield 
of  dried  algae  cake.  This  is  reasonable  because  the  more  dense  the  input 
solution,  the  faster  the  algae  deposits,  and  the  faster  the  web  velocity,  the 
greater  the  production  rate  of  dewatered  algae.  If  we  look  more  closely,  an 
increase  in  heating  energy  is  also  expected  as  the  web  velocity  increases  or  the 
concentration  of  input  solution  decreases,  given  the  desired  to  have  dried  algae 
cake  at  the  completion  of  the  HDD  process.  Table  5.7  demonstrates  the 
efficiency  of  the  overall  HDD  process.  As  can  be  observed,  an  increase  in  the 
concentration  of  the  input  solution  can  improve  the  efficiency  dramatically. 

According  to  Hooke's  Law,  and  Darcy's  Law,  changes  in  resistance  can  be 
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expected  when  there  exists  tension  in  the  web.  Moreover,  the  porosity,  pore 


sizes,  and  thickness  of  web  will  change.  The  fluid  dynamics  and  dryer 
requirements  are  therefore  going  to  be  different.  Hence,  an  alternative 
approach  to  HDD  control  is  necessary.  The  data  on  different  web  tensions,  input 
solution  concentrations,  and  web  velocities  is  needed  from  either  CFD 
simulations  or  experiments  with  the  process.  In  this  study,  CFD  modeling  and 
simulations  are  used  because  of  the  difficulty  in  getting  access  to  the  AVS  HDD 
process  and  the  appropriate  experimental  resources.  Tables  5.8~5.25  are  the 
relative  tables  estimated  by  CFD  simulations  for  different  web  tensions  and 
input  solution  concentrations.  As  the  tables  show,  the  energy  efficiency  can  be 
improved  especially  when  concentration  is  low  and  velocity  is  high.  The  possible 
reason  is  that  algae  do  not  dominate  the  resistance  in  those  conditions.  Fig.  5.5, 
shows  the  trend  of  energy  efficiency  with  different  input  solution 
concentrations  and  web  tension.  As  is  demonstrated,  the  higher  the  input 
solution  concentration,  the  higher  the  efficiency.  The  black  line  is  the  net  energy 
content  of  1  lb  of  dried  algae,  and  the  intersection,  at  about  12  g/L,  between 
the  concentration  curve  and  net  energy  content  curve  defines  the  operating 
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point  where  it  doesn't  make  sense  to  use  more  energy  to  make  products  with 


lower  energy  content.  Again,  if  we  compare  the  energy  efficiency  for  different 
web  tensions  at  10  g/L  concentration,  we  can  find  the  energy  efficiency  tradeoff 
point  at  15  N  web  tension.  If  we  compare  the  price  of  producing  dried  algae,  it 
follows  that  using  4.81  kBtu/hr  per  lb  of  dried  algae  is  the  best  energy  efficiency 
option  among  all  cases.  If  we  assume  that  the  algae  contains  60%  oil,  and  the 
density  of  oil  is  0.8  kg/m^  the  total  drying  energy  to  produce  1  gallon  of 
bio-diesel  will  be  about  19.59  kW.  If  the  energy  cost  is  8  cent  per  kW,  the  price 
to  produce  1  gallon  of  bio-diesel  will  be  1.56  US  dollar.  If  the  yield  of  algae  is 
5,000  gallons  per  year,  and  the  operating  cost  of  producing  algae  is  5,000  dollars 
per  acre  per  year,  the  selling  price  of  biodiesel  will  be  $2.89  plus  the  cost  to 
transport  the  algae,  extract  oil  from  the  algae,  and  profit.  The  cost  will  be  very 
competitive  if  the  concentration  of  algae  can  be  increased,  and  the  cost  of 
extracting  oil  from  the  algae  can  be  decreased.  Hence,  the  AVS  HDD  system 
should  be  an  efficient  way  to  generate  oil.  From  the  data  available  in  the  tables, 
appropriate  real-time  control  of  the  HDD  system  can  be  used  to  accommodate 
for  different  operating  situations  such  as  limitations  on  the  input  solution,  web 
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velocity,  heater  performance,  maximum  heat  supply,  and  so  on.  It  can  also  be 
used  to  adjust  the  operation  to  meet  specific  machine  or  operator  requirements, 
or  performance  restrictions  such  as  heater,  web  tension  limits,  pore  size  limits, 
or  material  limits.  Also,  the  tables  can  provide  important  data  to  estimate  the 
economic  cost  of  production  for  different  algae  yields,  including  the  cost  of 
preprocessing  the  input  solution  to  increase  concentration  or  the  modifying  the 
drying  method. 


Table  5.1  Expected  draining  Gravity  flow  rate  (mVs) 


weB'velogi^  (in/min) 
concentratioiTfe/Li 

1 

2 

4 

5 

6 

8 

9 

10 

5.8509E-06 

8.7766E-06 

1.1125E-05 

1.3164E-05 

1.4998E-05 

1.6685E-05 

1.8258E-05 

1.9739E-05 

2.1146E-05 

2.2488E-05 

2 

3.9002E-06 

5.8510E-06 

7.4174E-06 

8.7769E-06 

l.OOOlE-05 

1.1126E-05 

1.2175E-05 

1.3164E-05 

1.4102E-05 

1.4998E-05 

3.0763E-06 

4.6153E-06 

5.8511E-06 

6.9236E-06 

7.8892E-06 

8.7770E-06 

9.6051E-06 

1.0385E-05 

1.1126E-05 

1.1833E-05 

4 

2.5997E-06 

3.9002E-06 

4.9446E-06 

5.8512E-06 

6.6672E-06 

7.4176E-06 

8.1176E-06 

8.7770E-06 

9.4030E-06 

l.OOOlE-05 

5 

2.2814E-06 

3.4228E-06 

4.3394E-06 

5.1350E-06 

5.8512E-06 

6.5099E-06 

7.1242E-06 

7.7030E-06 

8.2525E-06 

8.7770E-06 

6 

2.0505E-06 

3.0764E-06 

3.9003E-06 

4.6154E-06 

5.2592E-06 

5.8512E-06 

6.4035E-06 

6.9238E-06 

7.4176E-06 

7.8892E-06 

1.8736E-06 

2.8110E-06 

3.5638E-06 

4.2173E-06 

4.8056E-06 

5.3466E-06 

5.8512E-06 

6.3267E-06 

6.7780E-06 

7.2089E-06 

1.7327E-06 

2.5997E-06 

3.2960E-06 

3.9003E-06 

4.4444E-06 

4.9448E-06 

5.4115E-06 

5.8512E-06 

6.2687E-06 

6.6672E-06 

9 

1.6174E-06 

2.4265E-06 

3.0764E-06 

3.6406E-06 

4.1484E-06 

4.6155E-06 

5.0511E-06 

5.4616E-06 

5.8512E-06 

6.2233E-06 

10 

1.5278E-06 

2.2814E-06 

2.8925E-06 

3.4229E-06 

3.9004E-06 

4.3395E-06 

4.7492E-06 

5.1351E-06 

5.5014E-06 

5.8512E-06 

11 

1.4382E-06 

2.1576E-06 

2.7355E-06 

3.2372E-06 

3.6888E-06 

4.1041E-06 

4.4915E-06 

4.8566E-06 

5.2031E-06 

5.5339E-06 

12 

1.3667E-06 

2.0505E-06 

2.5997E-06 

3.0765E-06 

3.5057E-06 

3.9004E-06 

4.2686E-06 

4.6155E-06 

4.9448E-06 

5.2592E-06 

13 

1.3042E-06 

1.9567E-06 

2.4808E-06 

2.9357E-06 

3.3453E-06 

3.7219E-06 

4.0733E-06 

4.4043E-06 

4.7185E-06 

5.0186E-06 

14 

1.2488E-06 

1.8736E-06 

2.3755E-06 

2.8111E-06 

3.2033E-06 

3.5640E-06 

3.9004E-06 

4.2174E-06 

4.5183E-06 

4.8056E-06 

15 

1.1994E-06 

1.7995E-06 

2.2815E-06 

2.6999E-06 

3.0765E-06 

3.4229E-06 

3.7461E-06 

4.0505E-06 

4.3395E-06 

4.6155E-06 

16 

1.1549E-06 

1.7328E-06 

2.1969E-06 

2.5998E-06 

2.9625E-06 

3.2961E-06 

3.6072E-06 

3.9004E-06 

4.1787E-06 

4.4444E-06 

17 

1.1146E-06 

1.6724E-06 

2.1203E-06 

2.5092E-06 

2.8592E-06 

3.1812E-06 

3.4815E-06 

3.7645E-06 

4.0331E-06 

4.2895E-06 

18 

1.0780E-06 

1.6173E-06 

2.0505E-06 

2.4266E-06 

2.7652E-06 

3.0765E-06 

3.3670E-06 

3.6406E-06 

3.9004E-06 

4.1484E-06 

19 

1.0444E-06 

1.5670E-06 

1.9867E-06 

2.3510E-06 

2.6790E-06 

2.9807E-06 

3.2621E-06 

3.5272E-06 

3.7789E-06 

4.0192E-06 

20 

1.0135E-06 

1.5206E-06 

1.9279E-06 

2.2815E-06 

2.5998E-06 

2.8926E-06 

3.1656E-06 

3.4229E-06 

3.6672E-06 

3.9004E-06 
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Table  5.2  Expected  draining  Capillary  flow  rate  (mVs) 


(in/min) 

concentrationTs/Li 

1 

2 

4 

5 

6 

7 

8 

9 

10 

1.3217E-06 

1.9833E-06 

2.5148E-06 

2.9762E-06 

3.3917E-06 

3.7739E-06 

4.1305E-06 

4.4664E-06 

4.7854E-06 

5.0900E-06 

2 

8.8081E-07 

1.3217E-06 

1.6758E-06 

1.9832E-06 

2.2600E-06 

2.5146E-06 

2.7522E-06 

2.9760E-06 

3.1885E-06 

3.3913E-06 

6.9468E-07 

1.0424E-06 

1.3216E-06 

1.5641E-06 

1.7824E-06 

1.9832E-06 

2.1705E-06 

2.3469E-06 

2.5145E-06 

2.6745E-06 

4 

5.8700E-07 

8.8079E-07 

1.1168E-06 

1.3216E-06 

1.5060E-06 

1.6757E-06 

1.8340E-06 

1.9832E-06 

2.1248E-06 

2.2600E-06 

5 

5.1512E-07 

7.7292E-07 

9.7999E-07 

1.1597E-06 

1.3216E-06 

1.4705E-06 

1.6094E-06 

1.7403E-06 

1.8645E-06 

1.9832E-06 

6 

4.6297E-07 

6.9467E-07 

8.8077E-07 

1.0423E-06 

1.1878E-06 

1.3216E-06 

1.4464E-06 

1.5641E-06 

1.6757E-06 

1.7824E-06 

4.2302E-07 

6.3472E-07 

8.0476E-07 

9.5237E-07 

1.0853E-06 

1.2075E-06 

1.3216E-06 

1.4291E-06 

1.5311E-06 

1.6285E-06 

3.9121E-07 

5.8699E-07 

7.4424E-07 

8.8075E-07 

1.0037E-06 

1.1167E-06 

1.2222E-06 

1.3216E-06 

1.4159E-06 

1.5060E-06 

9 

3.6515E-07 

5.4788E-07 

6.9465E-07 

8.2206E-07 

9.3677E-07 

1.0423E-06 

1.1408E-06 

1.2335E-06 

1.3216E-06 

1.4057E-06 

10 

3.4331E-07 

5.1510E-07 

6.5310E-07 

7.7289E-07 

8.8073E-07 

9.7995E-07 

1.0725E-06 

1.1597E-06 

1.2425E-06 

1.3216E-06 

11 

3.2594E-07 

4.8715E-07 

6.1765E-07 

7.3094E-07 

8.3293E-07 

9.2676E-07 

1.0143E-06 

1.0968E-06 

1.1751E-06 

1.2499E-06 

12 

3.0856E-07 

4.6296E-07 

5.8698E-07 

6.9464E-07 

7.9156E-07 

8.8073E-07 

9.6392E-07 

1.0423E-06 

1.1167E-06 

1.1878E-06 

13 

2.9443E-07 

4.4176E-07 

5.601 lE-07 

6.6284E-07 

7.5532E-07 

8.4041E-07 

9.1979E-07 

9.9458E-07 

1.0656E-06 

1.1334E-06 

14 

2.8193E-07 

4.2301E-07 

5.3632E-07 

6.3469E-07 

7.2325E-07 

8.0473E-07 

8.8073E-07 

9.5235E-07 

1.0203E-06 

1.0853E-06 

15 

2.7077E-07 

4.0626E-07 

5.1509E-07 

6.0957E-07 

6.9462E-07 

7.7287E-07 

8.4586E-07 

9.1465E-07 

9.7995E-07 

1.0423E-06 

16 

2.6073E-07 

3.9120E-07 

4.9599E-07 

5.8696E-07 

6.6886E-07 

7.4421E-07 

8.1450E-07 

8.8073E-07 

9.4361E-07 

1.0037E-06 

17 

2.5164E-07 

3.7756E-07 

4.7870E-07 

5.6650E-07 

6.4554E-07 

7.1826E-07 

7.8610E-07 

8.5002E-07 

9.1070E-07 

9.6865E-07 

18 

2.4336E-07 

3.6513E-07 

4.6295E-07 

5.4785E-07 

6.2430E-07 

6.9462E-07 

7.6022E-07 

8.2204E-07 

8.8073E-07 

9.3677E-07 

19 

2.3578E-07 

3.5376E-07 

4.4852E-07 

5.3078E-07 

6.0484E-07 

6.7298E-07 

7.3654E-07 

7.9643E-07 

8.5329E-07 

9.0758E-07 

20 

2.2880E-07 

3.4329E-07 

4.3525E-07 

5.1508E-07 

5.8695E-07 

6.5307E-07 

7.1475E-07 

7.7287E-07 

8.2804E-07 

8.8073E-07 

Table  5.3  Expected  output  algae  mass  flow  rate  (kg/s) 


weD^eleci^  (in/min) 
concentratioiTig4J 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

5.85673E-06 

8.78539E-06 

1.11366E-05 

1.31769E-05 

1.50133E-05 

1.67018E-05 

1.83E-05 

1.98E-05 

2.12E-05 

2.25E-05 

2 

7.81593E-06 

1.17254E-05 

1.48645E-05 

1.7589E-05 

2.00413E-05 

2.22966E-05 

2.44E-05 

2.64E-05 

2.83E-05 

3.01E-05 

3 

9.25677E-06 

1.38875E-05 

1.7606E-05 

2.08334E-05 

2.37388E-05 

2.64107E-05 

2.89E-05 

3.13E-05 

3.35E-05 

3.56E-05 

4 

1.04404E-05 

1.56635E-05 

1.9858E-05 

2.34986E-05 

2.6776E-05 

2.97902E-05 

3.26E-05 

3.53E-05 

3.78E-05 

4.02E-05 

5 

1.14642E-05 

1.71999E-05 

2.1806E-05 

2.58041E-05 

2.94033E-05 

3.27135E-05 

3.58E-05 

3.87E-05 

4.15E-05 

4.41E-05 

6 

1.23772E-05 

1.85697E-05 

2.35429E-05 

2.78596E-05 

3.17458E-05 

3.532E-05 

3.87E-05 

4.18E-05 

4.48E-05 

4.76E-05 

7 

1.32075E-05 

1.98157E-05 

2.51228E-05 

2.97294E-05 

3.38765E-05 

3.76908E-05 

4.12E-05 

4.46E-05 

4.78E-05 

5.08E-05 

8 

1.39736E-05 

2.09653E-05 

2.65804E-05 

3.14543E-05 

3.58422E-05 

3.9878E-05 

4.36E-05 

4.72E-05 

5.06E-05 

5.38E-05 

9 

1.46888E-05 

2.20371E-05 

2.79393E-05 

3.30627E-05 

3.7675E-05 

4.19173E-05 

4.59E-05 

4.96E-05 

5.31E-05 

5.65E-05 

10 

1.16371E-05 

2.30446E-05 

2.92169E-05 

3.45746E-05 

3.9398E-05 

4.38344E-05 

4.8E-05 

5.19E-05 

5.56E-05 

5.91E-05 

11 

1.59962E-05 

2.39981E-05 

3.04257E-05 

3.60052E-05 

4.10283E-05 

4.56459E-05 

5E-05 

5.4E-05 

5.79E-05 

6.16E-05 

12 

1.65996E-05 

2.4905E-05 

3.15757E-05 

3.73662E-05 

4.25792E-05 

4.7374E-05 

5.18E-05 

5.61E-05 

6.01E-05 

6.39E-05 

13 

1.71779E-05 

2.57717E-05 

3.26747E-05 

3.86667E-05 

4.40613E-05 

4.90231E-05 

5.37E-05 

5.8E-05 

6.22E-05 

6.61E-05 

14 

1.77314E-05 

2.66031E-05 

3.37287E-05 

3.99141E-05 

4.54828E-05 

5.06047E-05 

5.54E-05 

5.99E-05 

6.42E-05 

6.82E-05 

15 

1.8265E-05 

2.7403E-05 

3.4743E-05 

4.11145E-05 

4.68507E-05 

5.21268E-05 

5.7E-05 

6.17E-05 

6.61E-05 

7.03E-05 

16 

1.87789E-05 

2.8175E-05 

3.57218E-05 

4.22728E-05 

4.81707E-05 

5.35955E-05 

5.87E-05 

6.34E-05 

6.8E-05 

7.23E-05 

17 

1.92759E-05 

2.89217E-05 

3.66686E-05 

4.33933E-05 

4.94476E-05 

5.50162E-05 

6.02E-05 

6.51E-05 

6.98E-05 

7.42E-05 

18 

1.97597E-05 

2.96456E-05 

3.75864E-05 

4.44795E-05 

5.06854E-05 

5.63935E-05 

6.17E-05 

6.67E-05 

7.15E-05 

7.6E-05 

19 

2.02279E-05 

3.03487E-05 

3.84779E-05 

4.55346E-05 

5.18877E-05 

5.77313E-05 

6.32E-05 

6.83E-05 

7.32E-05 

7.79E-05 

20 

2.06837E-05 

3.10328E-05 

3.93453E-05 

4.6561  lE-05 

5.30575E-05 

5.90329E-05 

6.46E-05 

6.99E-05 

7.48E-05 

7.96E-05 
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Table  5.4  Expected  water  flow  rate  (kg/s)  in  the  output  of  capillary  zone 


(in/min) 

concentrationTs/Li 

1 

2 

4 

5 

6 

7 

8 

9 

10 

0.000502665 

0.000753322 

0.000954237 

0.00112836 

0.001284908 

0.001428635 

0.001562531 

0.001688543 

0.001808033 

0.001921997 

2 

0.000332614 

0.00049869 

0.000631913 

0.000747446 

0.000851379 

0.000946844 

0.001035834 

0.001119618 

0.001199099 

0.001274937 

0.000260308 

0.000390351 

0.000494706 

0.000585226 

0.000666681 

0.000741514 

0.00081128 

0.000876988 

0.000939328 

0.000998821 

4 

0.000218211 

0.000327258 

0.00041478 

0.000490711 

0.000559048 

0.000621803 

0.000680305 

0.000735395 

0.000787664 

0.000837545 

5 

0.000189935 

0.000284872 

0.000361079 

0.0004272 

0.000486712 

0.00054137 

0.000592328 

0.000640316 

0.000685849 

0.000729306 

6 

0.0001693 

0.000253936 

0.000321879 

0.000380838 

0.000433905 

0.000482644 

0.000528091 

0.000570889 

0.000611501 

0.000650262 

0.0001534 

0.000230095 

0.000291669 

0.000345103 

0.0003932 

0.000437378 

0.000478569 

0.000517367 

0.000554181 

0.000589319 

0.000140665 

0.000211001 

0.000267471 

0.000316479 

0.000360595 

0.000401115 

0.000438899 

0.000474485 

0.000508255 

0.000540492 

9 

0.000130189 

0.000195262 

0.000247525 

0.000292884 

0.000333715 

0.000371221 

0.000406193 

0.000439134 

0.000470393 

0.000500231 

10 

0.000123518 

0.000181996 

0.000230713 

0.000272995 

0.000311058 

0.000346021 

0.000378624 

0.000409333 

0.000438475 

0.000466293 

11 

0.000112505 

0.000170615 

0.000216289 

0.00025593 

0.000291618 

0.000324399 

0.000354968 

0.000383761 

0.000411086 

0.00043717 

12 

0.000107122 

0.000160708 

0.000203732 

0.000241075 

0.000274694 

0.000305575 

0.000334373 

0.000361498 

0.000387241 

0.000411814 

13 

0.000101323 

0.000151979 

0.000192668 

0.000227985 

0.000259781 

0.000288988 

0.000316225 

0.000341881 

0.000366229 

0.000389471 

14 

9.61296E-05 

0.000144209 

0.00018282 

0.000216333 

0.000246505 

0.000274223 

0.00030007 

0.000324417 

0.000347523 

0.00036958 

15 

9.14853E-05 

0.000137232 

0.000173976 

0.00020587 

0.000234585 

0.000260964 

0.000285562 

0.000308733 

0.000330724 

0.000351716 

16 

8.72658E-05 

0.000130919 

0.000165974 

0.000196403 

0.000223799 

0.000248966 

0.000272435 

0.000294542 

0.000315524 

0.000335552 

17 

8.34203E-05 

0.00012517 

0.000158687 

0.00018778 

0.000213975 

0.000238037 

0.000260478 

0.000281616 

0.000301678 

0.000320829 

18 

7.99302E-05 

0.000119902 

0.00015201 

0.000179881 

0.000204974 

0.000228025 

0.000249523 

0.000269773 

0.000288992 

0.000307339 

19 

7.66883E-05 

0.000115051 

0.000145861 

0.000172604 

0.000196685 

0.000218804 

0.000239434 

0.000258866 

0.000277309 

0.000294915 

20 

7.36987E-05 

0.000110563 

0.000140172 

0.000165873 

0.000189014 

0.000210273 

0.000230098 

0.000248774 

0.000266499 

0.000283419 

Inlet  air : 
Volume 
Humidity 
Temp 


Air  pre 
heater 


Heat 

supply 


Wet  algae 
Mass  flow  rate 
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content 
Temp 


Exhausted  gas 


Dryer 


Dried  algae 
cake 


Fig.  5.4  Schematic  for  algae  drying  process 
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Table  5.5  Required  airflow  for  heater  to  dry  algae  (kib/hr) 


(in/min) 

concentratioiTfs/Li 

1 

2 

4 

5 

6 

7 

8 

9 

10 

3.5580E-02 

5.3322E-02 

6.7543E-02 

7.9868E-02 

9.0949E-02 

1.0112E-01 

1.1060E-01 

1.1952E-01 

1.2798E-01 

1.3604E-01 

2 

2.3538E-02 

3.5290E-02 

4.4717E-02 

5.2893E-02 

6.0248E-02 

6.7004E-02 

7.3301E-02 

7.9230E-02 

8.4854E-02 

9.0221E-02 

1.8416E-02 

2.7617E-02 

3.4999E-02 

4.1403E-02 

4.7166E-02 

5.2460E-02 

5.7396E-02 

6.2045E-02 

6.6455E-02 

7.0664E-02 

4 

1.5434E-02 

2.3147E-02 

2.9337E-02 

3.4708E-02 

3.9542E-02 

4.3980E-02 

4.8118E-02 

5.2014E-02 

5.571  lE-02 

5.9240E-02 

5 

1.3431E-02 

2.0144E-02 

2.5533E-02 

3.0208E-02 

3.4416E-02 

3.8281E-02 

4.1885E-02 

4.5278E-02 

4.8498E-02 

5.1571E-02 

6 

1.1968E-02 

1.7952E-02 

2.2755E-02 

2.6923E-02 

3.0674E-02 

3.4120E-02 

3.7333E-02 

4.0358E-02 

4.3229E-02 

4.5969E-02 

1.0841E-02 

1.6262E-02 

2.0614E-02 

2.4390E-02 

2.7789E-02 

3.0912E-02 

3.3823E-02 

3.6565E-02 

3.9167E-02 

4.1650E-02 

9.9388E-03 

1.4908E-02 

1.8898E-02 

2.2361E-02 

2.5478E-02 

2.8341E-02 

3.1011E-02 

3.3525E-02 

3.5911E-02 

3.8189E-02 

9 

9.1961E-03 

1.3793E-02 

1.7484E-02 

2.0688E-02 

2.3572E-02 

2.6222E-02 

2.8692E-02 

3.1019E-02 

3.3227E-02 

3.5334E-02 

10 

8.7227E-03 

1.2852E-02 

1.6292E-02 

1.9278E-02 

2.1966E-02 

2.4435E-02 

2.6737E-02 

2.8906E-02 

3.0963E-02 

3.2928E-02 

11 

7.9422E-03 

1.2045E-02 

1.5269E-02 

1.8068E-02 

2.0587E-02 

2.2901E-02 

2.5059E-02 

2.7092E-02 

2.9021E-02 

3.0862E-02 

12 

7.5601E-03 

1.1342E-02 

1.4378E-02 

1.7014E-02 

1.9387E-02 

2.1566E-02 

2.3598E-02 

2.5513E-02 

2.7330E-02 

2.9064E-02 

13 

7.1487E-03 

1.0723E-02 

1.3594E-02 

1.6085E-02 

1.8329E-02 

2.0389E-02 

2.2311E-02 

2.4121E-02 

2.5839E-02 

2.7479E-02 

14 

6.7803E-03 

1.0171E-02 

1.2895E-02 

1.5259E-02 

1.7387E-02 

1.9342E-02 

2.1165E-02 

2.2882E-02 

2.4512E-02 

2.6067E-02 

15 

6.4507E-03 

9.6763E-03 

1.2267E-02 

1.4516E-02 

1.6541E-02 

1.8401E-02 

2.0135E-02 

2.1769E-02 

2.3320E-02 

2.4800E-02 

16 

6.1512E-03 

9.2283E-03 

1.1699E-02 

1.3844E-02 

1.5775E-02 

1.7549E-02 

1.9203E-02 

2.0762E-02 

2.2241E-02 

2.3652E-02 

17 

5.8782E-03 

8.8201E-03 

1.1182E-02 

1.3232E-02 

1.5078E-02 

1.6773E-02 

1.8355E-02 

1.9844E-02 

2.1258E-02 

2.2607E-02 

18 

5.6304E-03 

8.4461E-03 

1.0708E-02 

1.2671E-02 

1.4439E-02 

1.6063E-02 

1.7577E-02 

1.9003E-02 

2.0357E-02 

2.1649E-02 

19 

5.4002E-03 

8.1017E-03 

1.0271E-02 

1.2154E-02 

1.3850E-02 

1.5408E-02 

1.6860E-02 

1.8229E-02 

1.9528E-02 

2.0767E-02 

20 

5.1879E-03 

7.7829E-03 

9.8672E-03 

1.1676E-02 

1.3305E-02 

1.4802E-02 

1.6197E-02 

1.7512E-02 

1.8760E-02 

1.9951E-02 

Table  5.6  Heating  energy  requirements  for  heating  air  (kBtu/hr) 


(in/min) 

concentratioiTts/Li 

1 

2 

4 

5 

6 

7 

8 

9 

10 

5.7659E+00 

8.6410E+00 

1.0946E+01 

1.2943E+01 

1.4739E+01 

1.6387E+01 

1.7923E+01 

1.9369E+01 

2.0739E+01 

2.2046E+01 

2 

3.8150E+00 

5.7199E+O0 

7.2479E+00 

8.5731E+00 

9.7652E+O0 

1.0860E+01 

1.1881E+01 

1.2842E+01 

1.3753E+01 

1.4623E+01 

2.9855E+00 

4.4770E+00 

5.6738E+00 

6.7120E+00 

7.6462E+O0 

8.5045E-f00 

9.3046E+00 

1.0058E+01 

1.0773E+01 

1.1456E+01 

4 

2.5025E+00 

3.7531E+O0 

4.7568E+00 

5.6276E+00 

6.4113E+O0 

7.1310E+00 

7.8020E+00 

8.4337E+00 

9.0332E+00 

9.6052EH-00 

5 

2.1781E+00 

3.2668E+O0 

4.1407E+00 

4.8989E+00 

5.5814E+O0 

6.2082E+00 

6.7925E+00 

7.3428E+00 

7.8650E+00 

8.3633E+00 

6 

1.9413E+00 

2.9118E+O0 

3.6909E+00 

4.3670E+00 

4.9755E+O0 

5.5344E4-00 

6.0555E+00 

6.5462E+00 

7.0119E+00 

7.4564E+00 

1.7589E+00 

2.6383E+00 

3.3443E+00 

3.9569E+00 

4.5084E+00 

5.0149E+00 

5.4872E+00 

5.9321E+00 

6.3542E+00 

6.7571E+00 

1.6127E+00 

2.4191E+O0 

3.0666E+00 

3.6285E+00 

4.1343E+O0 

4.5988E4-00 

5.0320E+00 

5.4400E+00 

5.8272E+00 

6.1968E+00 

9 

1.4925E+00 

2.2385E+O0 

2.8377E+00 

3.3577E+00 

3.8258E+O0 

4.2558E-f00 

4.6567E+00 

5.0343E+00 

5.3927E+00 

5.7348E+00 

10 

1.4159E+00 

2.0863E+00 

2.6448E+00 

3.1294E+00 

3.5658E+O0 

3.9666E4-00 

4.3403E+00 

4.6923E+00 

5.0264E+00 

5.3453E+00 

11 

1.2896E+00 

1.9557E+00 

2.4792E+00 

2.9336E+00 

3.3427E+00 

3.7184E+00 

4.0688E+00 

4.3989E+00 

4.7121E+00 

5.0111E+00 

12 

1.2278E+00 

1.8420E+00 

2.3351E+00 

2.7631E+00 

3.1484E+00 

3.5024E+00 

3.8325E+00 

4.1434E+00 

4.4384E+00 

4.7201E+00 

13 

1.1612E+00 

1.7418E+00 

2.2081E+00 

2.6129E+00 

2.9773E+00 

3.3120E+00 

3.6242E+00 

3.9182E+00 

4.1972E+00 

4.4636E+00 

14 

1.1016E+00 

1.6526E+00 

2.095  lE+00 

2.4791E+00 

2.8249E+00 

3.1425E+00 

3.4387E+00 

3.7177E+00 

3.9825E+00 

4.2353E+00 

15 

1.0483E+00 

1.5725E+00 

1.9936E+00 

2.3590E+00 

2.6881E+00 

2.9903E+00 

3.2722E+00 

3.5377E+00 

3.7897E+00 

4.0303E+00 

16 

9.9988E-01 

1.5001E+00 

1.9017E+00 

2.2504E+00 

2.5643E+00 

2.8526E+00 

3.1215E+00 

3.3748E+00 

3.6152E+00 

3.8447E+00 

17 

9.5573E-01 

1.4341E+00 

1.8180E+00 

2.1514E+00 

2.4515E+00 

2.7272E+00 

2.9843E+00 

3.2264E+00 

3.4563E+00 

3.6757E+00 

18 

9.1567E-01 

1.3736E+00 

1.7414E+00 

2.0607E+00 

2.3482E+00 

2.6122E+00 

2.8585E+00 

3.0905E+00 

3.3106E+00 

3.5208E+00 

19 

8.7845E-01 

1.3179E+00 

1.6708E+00 

1.9772E+00 

2.2530E+00 

2.5064E+00 

2.7427E+00 

2.9653E+00 

3.1765E+00 

3.3782E+00 

20 

8.4413E-01 

1.2664E+00 

1.6055E+00 

1.8999E+00 

2.1649E+00 

2.4084E+00 

2.6355E+00 

2.8494E+00 

3.0524E+00 

3.2462E+00 
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Table  5.7  Energy  efficiency  (kBtu/lb  of  dried  algae) 


(in/min) 

concentrationTs/Li 

1 

2 

4 

5 

6 

7 

8 

9 

10 

1.2415E+05 

1.2404E+05 

1.2395E+05 

1.2387E+05 

1.2380E+05 

1.2374E+05 

1.2368E+05 

1.2362E+05 

1.2356E+05 

1.2351E+05 

2 

6.1556E+04 

6.1520E+04 

6.1492E+04 

6.1468E+04 

6.1448E+04 

6.1427E+04 

6.1408E+04 

6.1390E+04 

6.1373E+04 

6.1356E+04 

4.0673E+04 

4.0655E+04 

4.0641E+04 

4.0630E+04 

4.0620E+04 

4.0610E+04 

4.0600E+04 

4.0591E+04 

4.0583E+04 

4.0575E+04 

4 

3.0228E+04 

3.0217E+04 

3.0209E+04 

3.0202E+04 

3.0197E+04 

3.0188E+04 

3.0181E+04 

3.0174E+04 

3.0167E+04 

3.0160E+04 

5 

2.3960E+04 

2.3952E+04 

2.3947E+04 

2.3942E+04 

2.3939E+04 

2.3933E4-04 

2.3928E+04 

2.3923E+04 

2.3918E+04 

2.3913E+04 

6 

1.9780E+04 

1.9775E+04 

1.9771E+04 

1.9768E+04 

1.9765E+04 

1.9761E4-04 

1.9757E+04 

1.9753E+04 

1.9750E+04 

1.9746E+04 

1.6794E+04 

1.6790E+04 

1.6787E+04 

1.6785E+04 

1.6783E+04 

1.6780E4-04 

1.6777E+04 

1.6774E+04 

1.6771E+04 

1.6768E-h04 

1.4555E+04 

1.4552E+04 

1.4549E+04 

1.4548E+04 

1.4546E+04 

1.4544E4-04 

1.4541E+04 

1.4539E+04 

1.4536E+04 

1.4534E+04 

9 

1.2814E+04 

1.2810E+04 

1.2809E+04 

1.2807E+04 

1.2806E+04 

1.2804E4-04 

1.2802E+04 

1.2800E+04 

1.2798E+04 

1.2796E+04 

10 

1.1571E+04 

1.1417E+04 

1.1416E+04 

1.1415E+04 

1.1414E+04 

1.1412E4-04 

1.1410E+04 

1.1409E+04 

1.1407E+04 

1.1405E+04 

11 

1.0167E+04 

1.0277E+04 

1.0276E+04 

1.0275E+04 

1.0275E+04 

1.0273E4-04 

1.0271E+04 

1.0270E+04 

1.0269E+04 

1.0267E+04 

12 

9.3278E+03 

9.3272E+03 

9.3262E+03 

9.3255E+03 

9.3250E+03 

9.3236E+03 

9.3223E+03 

9.3210E+03 

9.3198E+03 

9.3187E+03 

13 

8.5251E+03 

8.5232E+03 

8.5224E+03 

8.5218E+03 

8.5215E+03 

8.5202E+03 

8.5191E+03 

8.5180E+03 

8.5169E+03 

8.5159E+03 

14 

7.8351E+03 

7.8341E+03 

7.8335E+03 

7.8330E+03 

7.8326E+03 

7.8316E+03 

7.8305E+03 

7.8296E+03 

7.8286E+03 

7.8277E+03 

15 

7.2381E+03 

7.2369E+03 

7.2363E+03 

7.2359E+03 

7.2356E+03 

7.2347E+03 

7.2338E+03 

7.2329E+03 

7.2320E+03 

7.2312E+03 

16 

6.7148E+03 

6.7142E+03 

6.7137E+03 

6.7134E+03 

6.7132E+03 

6.7124E+03 

6.7115E+03 

6.7107E+03 

6.7100E+03 

6.7093E+03 

17 

6.2528E+03 

6.2531E+03 

6.2527E+03 

6.2524E+03 

6.2522E+03 

6.2514E+03 

6.2507E+03 

6.2500E+03 

6.2493E+03 

6.2486E+03 

18 

5.8440E+03 

5.8432E+03 

5.8428E+03 

5.8426E+03 

5.8425E+03 

5.8417E+03 

5.841  lE+03 

5.8404E+03 

5.8398E+03 

5.8392E+03 

19 

5.4767E+03 

5.4764E+03 

5.4761E+03 

5.4759E+03 

5.4758E+03 

5.4751E+03 

5.4745E+03 

5.4739E+03 

5.4733E+03 

5.4728E+03 

20 

5.1468E+03 

5.1462E+03 

5.1460E+03 

5.1458E+03 

5.1458E+03 

5.1451E+03 

5.1446E+03 

5.1440E+03 

5.1435E+03 

5.1430E+03 

Table  5.8  Expected  draining  Gravity  flow  rate  (m^3/s)  for  different  web  tension  (1 

g/L  input  solution) 


""'■'■•Aaieb  velocity 
( i  n  /  rn  i  1  a  t  e  d 

height''''^.. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

5.8509E-06  8.7766E-06 

5.9094E-06  8.8644E-06 

1.1125E-05 

1.1237E-05 

1.3164E-05 

1.3295E-05 

1.4998E-05 

1.5148E-05 

1.6685E-05 

1.6852E-05 

1.8258E-05 

1.8440E-05 

1.9739E-05 

1.9937E-05 

2.1146E-05 

2.1357E-05 

2.2488E-05 

2.2713E-05 

5 

10 

5.9386E-06  8.9083E-06 

1.1292E-05 

1.3361E-05 

1.5223E-05 

1.6935E-05 

1.8532E-05 

2.0035E-05 

2.1463E-05 

2.2825E-05 

15 

5.9562E-06  8.9346E-06 

1.1326E-05 

1.3401E-05 

1.5268E-05 

1.6985E-05 

1.8586E-05 

2.0095E-05 

2.1526E-05 

2.2893E-05 

20 

5.9679E-06  8.9521E-06 

1.1348E-05 

1.3427E-05 

1.5298E-05 

1.7019E-05 

1.8623E-05 

2.0134E-05 

2.1569E-05 

2.2938E-05 

Table  5.9  Expected  draining  Capillary  flow  rate  (m^3/s)  for  different  web  tension  (1 

g/L  input  solution) 


velocity 

(in/niiny'''''''sU^^ 

height^'^...^ 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

1.3217E-06  1.9833E-06 

1.3307E-06  2.0009E-06 

2.5148E-06 

2.5399E-06 

2.9762E-06 

3.0081E-06 

3.3917E-06 

3.4300E-06 

3.7739E-06 

3.S181E-06 

4.1305E-06 

4.1803E-06 

4.4664E-06 

4.5217E-06 

4.7854E-06 

4.8459E-06 

5.0900E-06 

5.1555E-06 

5 

10 

1.3425E-06  2.0219E-06 

2.5689E-06 

3.0445E-06 

3.4730E-06 

3.8674E-06 

4.2355E-06 

4.5826E-06 

4.9122E-06 

5.2270E-06 

15 

1.3554E-06  2.0444E-06 

2.5997E-06 

3.0827E-06 

3.5181E-06 

3.9190E-06 

4.2932E-06 

4.6461E-06 

4.9S13E-06 

5.3014E-06 

20 

1.3688E-06  2.0677E-06 

2.6315E-06 

3.1221E-06 

3.5646E-06 

3.9720E-06 

4.3525E-06 

4.7113E-06 

5.0521E-06 

5.3777E-06 
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Table  5.10  Expected  output  algae  mass  flow  rate  (kg/s)  for  different  web  tension  (1 

g/L  input  solution) 


velocity 

( i  n  /  IT1  i  1  a  t  e  d 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

5.8567E-06  8.7854E-06 

1.1137E-05 

1.3177E-05 

1.5013E-05 

1.6702E-05 

1.8276E-05 

1.9759E-05 

2.1167E-05 

2.2510E-05 

5 

5.9153E-06  8.8732E-06 

1.1248E-05 

1.3309E-05 

1.5163E-05 

1.6869E-05 

1.8459E-05 

1.9957E-05 

2.1378E-05 

2.2736E-05 

10 

5.9446E-06  8.9172E-06 

1.1304E-05 

1.3375E-05 

1.5238E-05 

1.6952E-05 

1.8550E-05 

2.0056E-05 

2.1484E-05 

2.2848E-05 

15 

5.9622E-06  8.9435E-06 

1.1337E-05 

1.3414E-05 

1.5283E-05 

1.7002E-05 

1.8605E-05 

2.0115E-05 

2.1548E-05 

2.2916E-05 

20 

5.9739E-06  8.9611E-06 

1.1359E-05 

1.3440E-05 

1.5314E-05 

1.7036E-05 

1.8642E-05 

2.0154E-05 

2.1590E-05 

2.2961E-05 

Table  5.11  Airflow  (kib/hr)  requirements  for  heater  to  dry  algae  for  different  web 

tension  (1  g/L  input  solution) 


velocity 

(in^fujT) 

concentrationfe/O 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

0.0356  0.0533 

0.0675 

0.0799 

0.0909 

0.1011 

0.1106 

0.1195 

0.1280 

0.1360 

5 

0.0362  0.0540 

0.0682 

0.0805 

0.0916 

0.1017 

0.1111 

0.1200 

0.1284 

0.1364 

10 

0.0360  0.0535 

0.0674 

0.0794 

0.0902 

0.1000 

0.1092 

0.1178 

0.1260 

0.1338 

15 

0.0355  0.0525 

0.0660 

0.0776 

0.0880 

0.0975 

0.1063 

0.1146 

0.1225 

0.1300 

0.1256 

20 

0.0348  0.0512 

0.0642 

0.0754 

0.0853 

0.0945 

0.1029 

0.1109 

0.1184 

Table  5.12  Heating  energy  requirements  for  heating  air  (kBtu/hr)  for  different  web 

tension  (1  g/L  input  solution) 


velocity 

(iriMin) 

CO  n  ce  n  t  ra  1 1  o  rHg/Q 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

5.7659  8.6410 

10.9457 

12.9430 

14.7386 

16.3873 

17.9231 

19.3686 

20.7392 

22.0464 

5 

5.8714  8.7535 

11.0561 

13.0480 

14.8365 

16.4771 

18.0043 

19.4407 

20.8021 

22.1000 

10 

5.8412  8.6692 

10.9217 

12.8668 

14.6114 

16.2101 

17.6973 

19.0953 

20.4196 

21.6817 

15 

5.7566  8.5053 

10.6878 

12.5691 

14.2543 

15.7971 

17.2312 

18.5784 

19.8540 

21.0690 

20 

5.6445  8.3009 

10.4031 

12.2114 

13.8290 

15.3083 

16.6822 

17.9721 

19.1926 

20.3545 

Table  5.13  Energy  efficiency  (kBtu/lb  of  dried  algae)  for  different  web  tension  (1 

g/L  input  solution) 


^■''sveb  velocity 
(m/min) 

CO  n  ce  n  t  ra  t  i  0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

1.2415E+05 

1.2404E+05 

1.2441E+05 

1.2260E+05 

1.2395E+05 

1.2396E+05 

1.2185E+05 

1.2387E+05 

1.2364E+05 

1.2132E+05 

1.2380E+05  1.2374E+05 

1.2339E+05  1.2318E+05 

1.2092E+05  1.2059E+05 

1.2368E+05 

1.2301E+05 

1.2031E+05 

1.2362E+05  1.2356E+05 

1.2285E+05  1.2271E+05 

1.2007E+05  1.1986E+05 

1.2351E+05 

1.2259E+05 

1.1967E+05 

5 

1.2518E+05 

1.2392E+05 

10 

15 

1.2176E+05 

1.1993E+05 

1.1889E+05 

1.1817E+05 

1.1762E+05  1.1717E+05 

1.1680E+05 

1.164SE+05  1.1620E+05 

1.1595E+05 

20 

1.1916E+05 

1.1682E+05 

1.1549E+05 

1.1458E+05 

1.1389E+05  1.1332E+05 

1.1286E+05 

1.1246E+05  1.1211E+05 

1.1180E+05 
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Table  5.14  Expected  draining  Gravity  flow  rate  (mVs)  for  different  web  tension  (10 

g/L  input  solution) 


■^vveb  velocity 
( i  n/rTiiiT]''''''SMIuj|  ated 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

1.5278E-06 

2.2814E-06 

2.S925E-06 

3.4229E-06 

3.9004E-06  4.3395E-06 

4.7492E-06 

5.1351E-06  5.5014E-06 

5.8512E-06 

5 

1.5385E-06 

2.2974E-06 

2.9127E-06 

3.4468E-06 

3.9277E-06  4.3699E-06 

4.7824E-06 

5.1710E-06  5.5400E-06 

5.8922E-06 

10 

1.5477E-06 

2.3111E-06 

2.9301E-06 

3.4674E-06 

3.9511E-06  4.3959E-06 

4.8109E-06 

5.2018E-06  5.5730E-06 

5.9273E-06 

15 

1.5553E-06 

2.3225E-06 

2.9445E-06 

3.4845E-06 

3.9706E-06  4.4176E-06 

4.8346E-06 

5.2275E-06  5.6005E-06 

5.9566E-06 

20 

1.5614E-06 

2.3316E-06 

2.9561E-06 

3.4982E-06 

3.9862E-06  4.4350E-06 

4.8536E-06 

5.2481E-06  5.6225E-06 

5.9800E-06 

Table  5.15  Expected  draining  Capillary  flow  rate  (mVs)  for  different  web  tension 

(10  g/L  input  solution) 


velocity 

( i  n/rri  i  n)''''''sUiTU  1  ate^ 
height 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

3.4331E-07 

5.1510E-07 

6.5310E-07 

6.5626E-07 

7.7289E-07 

7.7724E-07 

8.8073E-07  9.7995E-07 

8.8623E-07  9.8655E-07 

1.0725E-06 

1.0802E-06 

1.1597E-06  1.2425E-06 

1.1684E-06  1.2523E-06 

1.3216E-06 

1.3323E-06 

5 

3.4398E-07 

5.1704E-07 

10 

3.4485E-07 

5.1923E-07 

6.5972E-07 

7.8191E-07 

8.9207E-07  9.9353E-07 

1.0883E-06 

1.1776E-06  1.2625E-06 

1.3435E-06 

15 

3.4590E-07 

5.2166E-07 

6.6346E-07 

7.8691E-07 

8.9829E-07  1.0009E-06 

1.0968E-06 

1.1S72E-06  1.2731E-06 

1.3552E-06 

20 

3.4715E-07 

5.2435E-07 

6.6752E-07 

7.9227E-07 

9.0488E-07  1.0087E-06 

1.1057E-06 

1.1973E-06  1.2843E-06 

1.3675E-06 

Table  5.16  Expected  output  algae  mass  flow  rate  (kg/s)  for  different  web  tension 

(10  g/L  input  solution) 


"'■■'■-^aieb  velocity 
( i  n  /  m  i  ny'''''5«3TU  1  a  t  e  d 
height 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

1.5432E-05  2.3045E-05 

2.9217E-05 

3.4575E-05 

3.9398E-05 

4.3834E-05 

4.7971E-05 

5.1870E-05 

5.5570E-05 

5.9104E-05 

5 

1.5540E-05  2.3206E-05 

2.9421E-05 

3.4817E-05 

3.9674E-05 

4.4140E-05 

4.8307E-05 

5.2233E-05 

5.5959E-05 

5.9517E-05 

10 

1.5633E-05  2.3344E-05 

2.9597E-05 

3.5024E-05 

3.9910E-05 

4.4404E-05 

4.8595E-05 

5.2544E-05 

5.6293E-05 

5.9872E-05 

15 

1.5710E-05  2.3459E-05 

2.9743E-05 

3.5197E-05 

4.0107E-05 

4.4623E-05 

4.8S35E-05 

5.2803E-05 

5.6570E-05 

6.0167E-05 

20 

1.5772E-05  2.3552E-05 

2.9860E-05 

3.5335E-05 

4.0265E-05 

4.4798E-05 

4.9027E-05 

5.3011E-05 

5.6793E-05 

6.0404E-05 

Table  5.17  Airflow  (kib/hr)  requirements  for  heater  to  dry  algae  for  different  web 

tension  (10  g/L  input  solution) 


''-■•sateb  velocity 
(mAwin) 

CO n ce n t ra t i 0 n  {^/U 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

8.7227E-03 

8.9061E-03 

1.2852E-02 

1.3060E-02 

1.6292E-02 

1.6506E-02 

1.9278E-02 

1.9488E-02 

2.1966E-02  2.4435E-02 

2.2167E-02  2.4624E-02 

2.6737E-02 

2.6912E-02 

2.8906E-02  3.0963E-02 

2.9065E-02  3.1106E-02 

3.2928E-02 

3.3053E-02 

5 

10 

9.0433E-03 

1.3201E-02 

1.6636E-02 

1.9601E-02 

2.2259E-02  2.4692E-02 

2.6956E-02 

2.9083E-02  3.1098E-02 

3.3019E-02 

15 

9.1339E-03 

1.3276E-02 

1.6684E-02 

1.9617E-02 

2.2241E-02  2.4639E-02 

2.6867E-02 

2.8959E-02  3.0938E-02 

3.2823E-02 

20 

9.1777E-03 

1.3282E-02 

1.6647E-02 

1.9534E-02 

2.2111E-02  2.4463E-02 

2.6644E-02 

2.8690E-02  3.0624E-02 

3.2463E-02 
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Table  5.18  Heating  energy  (kBtu/hr)  requirements  for  heating  air  for  different  web 

tension  (10  g/L  input  solution) 


velocity 

(m/min) 

CO  n  ce  n  t  rat  i  0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

1.4159 

2.0863 

2.6448 

3.1294 

3.5658  3.9666 

4.3403 

4.6923  5.0264 

5.3453 

5 

1.4457 

2.1200 

2.6794 

3.1635 

3.5984  3.9973 

4.3687 

4.7182  5.0496 

5.3656 

10 

1.4679 

2.1430 

2.7006 

3.1819 

3.6134  4.0084 

4.3759 

4.7212  5.0484 

5.3601 

15 

1.4826 

2.1550 

2.7083 

3.1844 

3.6105  3.9999 

4.3616 

4.7011  5.0225 

5.3284 

20 

1.4897 

2.1561 

2.7023 

3.1710 

3.5895  3.9713 

4.3255 

4.6576  4.9716 

5.2703 

Table  5.19  Energy  efficiency  (kBtu/lb  of  dried  algae)  for  different  web  tension  (10 

g/L  input  solution) 


velocity 

(rnAniti) 

CO n ce n t rat i 0 n  {&/U 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

1.1571E+04 

1.1732E+04 

1.1417E+04 

1.1521E+04 

1.1416E+04 

1.1485E+04 

1.1415E+04 

1.1459E+04 

1.1414E+04  1.1412E+04 

1.1438E+04  1.1420E+04 

1.1410E+04 

1.1405E+04 

1.1409E+04  1.1407E+04 

1.1392E+04  1.1380E+04 

1.1405E+04 

1.1369E+04 

5 

10 

1.1842E+04 

1.1577E+04 

1.1507E+04 

1.1457E+04 

1.1418E+04  1.1384E+04 

1.1356E+04 

1.1331E+04  1.1310E+04 

1.1290E+04 

15 

1.1902E+04 

1.1912E+04 

1.1585E+04 

1.1545E+04 

1.1483E+04 

1.1413E+04 

1.1410E+04 

1.1317E+04 

1.1353E+04  1.1304E+04 

1.1242E+04  1.1180E+04 

1.1263E+04 

1.1126E+04 

1.122SE+04  1.1196E+04 

1.1080E+04  1.1040E+04 

1.1168E+04 

1.1003E+04 

20 

Table  5.20  Expected  draining  Gravity  flow  rate  (mVs)  for  different  web  tension  (20 

g/L  input  solution) 


^■-'SAieb  velocity 
(in/rniiT)''''''''siQ^^ 

height^\..^ 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

1.0135E-06 

1.5206E-06  1.9279E-06 

2.2815E-06  2.5998E-06 

2.8926E-06 

3.1656E-06 

3.4229E-06 

3.6672E-06 

3.9004E-06 

5 

1.0135E-06 

1.5206E-06  1.9279E-06 

2.2815E-06  2.5998E-06 

2.8926E-06 

3.1656E-06 

3.1815E-06 

3.4229E-06 

3.4401E-06 

3.6672E-06 

3.6855E-06 

3.9004E-06 

3.9199E-06 

10 

1.0186E-06 

1.5282E-06  1.9376E-06 

2.2929E-06  2.6128E-06 

2.9070E-06 

15 

1.0226E-06 

1.5343E-06  1.9453E-06 

2.3020E-06  2.6232E-06 

2.9186E-06 

3.1941E-06 

3.4537E-06 

3.7002E-06 

3.9355E-06 

20 

1.0257E-06 

1.5389E-06  1.9511E-06 

2.3089E-06  2.6310E-06 

2.9273E-06 

3.2036E-06 

3.4640E-06 

3.7112E-06 

3.9472E-06 

Table  5.21  Expected  draining  Capillary  flow  rate  (mVs)  for  different  web  tension 

(20  g/L  input  solution) 


"-'^j/veb  velocity 
( i  ri/rn  i  n)‘'''-ski!ul  ate  d 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

2.2880E-07 

2.3004E-07 

3.4329E-07  4.3525E-07 

3.4556E-07  4.3846E-07 

5.1508E-07  5.8695E-07 

5.1917E-07  5.9186E-07 

6.5307E-07 

6.5877E-07 

7.1475E-07 

7.2121E-07 

7.7287E-07 

7.8006E-07 

8.2804E-07 

8.3594E-07 

8.8073E-07 

8.8932E-07 

5 

10 

2.3044E-07 

3.4667E-07  4.402SE-07 

5.2168E-07  5.9505E-07 

6.6261E-07 

7.2569E-07 

7.S517E-07 

8.4168E-07 

8.9566E-07 

15 

2.3098E-07 

3.4798E-07  4.4233E-07 

5.2445E-07  5.9852E-07 

6.6677E-07 

7.3051E-07 

7.9064E-07 

8.4777E-07 

9.0237E-07 

20 

2.3167E-07 

3.4948E-07  4.4462E-07 

5.2749E-07  6.0229E-07 

6.7124E-07 

7.3567E-07 

7.9647E-07 

8.5426E-07 

9.0949E-07 
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Table  5.22  Expected  output  algae  mass  flow  rate  (kg/s)  for  different  web  tension 

(20  g/L  input  solution) 


velocity 

( i  n/tTi  i  ate  d 

height 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

2.0684E-05 

3.1033E-05  3.9345E-05 

4.6561E-05  5.3058E-05 

5.9032E-05 

6.4605E-05 

6.9856E-05 

7.4841E-05 

7.9600E-05 

4 

2.0684E-05 

3.1033E-05  3.9345E-05 

4.6561E-05  5.3058E-05 

5.9032E-05 

6.4605E-05 

6.9856E-05 

7.4841E-05 

7.9600E-05 

8 

2.0787E-05 

3.1188E-05  3.9542E-05 

4.6794E-05  5.3323E-05 

5.9327E-05 

6.4928E-05 

7.0205E-05 

7.5215E-05 

7.9998E-05 

12 

2.0870E-05 

2.0932E-05 

3.1312E-05  3.9699E-05 

3.1405E-05  3.9817E-05 

4.6980E-05  5.3535E-05 

4.7120E-05  5.3694E-05 

5.9563E-05 

5.9740E-05 

6.5186E-05 

6.5380E-05 

7.0485E-05 

7.0694E-05 

7.5514E-05 

7.5739E-05 

8.0316E-05 

8.0555E-05 

16 

Table  5.23  Airflow  (kib/hr)  requirements  for  heater  to  dry  algae  for  different  web 

tension  (20  g/L  input  solution) 


"^ifi/eb  velocity 

CO n ce n trati 0 n 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

5.1879E-03 

7.7829E-03  9.8672E-03 

1.1676E-02  1.3305E-02 

1.4802E-02  1.6197E-02 

1.7512E-02 

1.8760E-02 

1.9951E-02 

5 

5.1000E-03 

7.6221E-03  9.6400E-03 

1.1387E-02  1.2958E-02 

1.4398E-02  1.5740E-02 

1.7003E-02 

1.8201E-02 

1.9343E-02 

10 

5.1789E-03 

7.7041E-03  9.7148E-03 

1.1450E-02  1.3007E-02 

1.4431E-02  1.5756E-02 

1.7002E-02 

1.8182E-02 

1.9306E-02 

15 

5.2261E-03 

7.7401E-03  9.7325E-03 

1.1447E-02  1.2980E-02 

1.4381E-02  1.5683E-02 

1.6904E-02 

1.8059E-02 

1.9159E-02 

20 

5.2415E-03 

7.7299E-03  9.6926E-03 

1.1376E-02  1.2878E-02 

1.4248E-02  1.5518E-02 

1.6708E-02 

1.7833E-02 

1.8902E-02 

Table  5.24  Heating  energy  (kBtu/hr)  requirements  for  heating  air  for  different  web 

tension  (20  g/L  input  solution) 


^■-'Sftieb  velocity 
(iriTmyi) 

conce  ntratiorifg/t) 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

0.8441 

1.2664  1.6055 

1.8999  2.1649 

2.4084 

2.6355 

2.8494 

3.0524 

3.2462 

5 

0.8299 

0.8427 

1.2403  1.5687 

1.2536  1.5808 

1.8530  2.1086 

1.8633  2.1165 

2.3430 

2.3484 

2.5614 

2.5641 

2.7669 

2.7668 

2.9618 

2.9588 

3.1477 

3.1417 

10 

15 

0.8503 

1.2595  1.5837 

1.8627  2.1123 

2.3404 

2.5522 

2.7510 

2.9390 

3.1181 

20 

0.8528 

1.2578  1.5773 

1.8513  2.0958 

2.3188 

2.5255 

2.7192 

2.9023 

3.0764 

Table  5.25  Energy  efficiency  (kBtu/lb  of  dried  algae)  for  different  web  tension  (20 

g/L  input  solution) 


velocity 

(inTmin) 

CO  n  ce  n  trati  0  n^g/jJ 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

5.1468E+03 

5.1462E+03  5.1460E+03 

5.1458E+03  5.1458E+03 

5.1451E+03  5.1446E+03 

5.1440E+03 

5.1435E+03 

5.1430E+03 

5 

5.0599E+03 

5.0403E+03  5.0280E+03 

5.0190E+03  5.0119E+03 

5.0055E+03  5.0000E+03 

4.9951E+03 

4.9908E+03 

4.9S69E+03 

10 

5.1124E+03 

5.0691E+03  5.0417E+03 

5.0216E+03  5.0057E+03 

4.9921E+03  4.9803E+03 

4.9701E+03 

4.9609E+03 

4.9527E+03 

15 

5.1384E+03 

5.0726E+03  5.0309E+03 

5.0002E+03  4.9759E+03 

4.9552E+03  4.9375E+03 

4.9220E+03 

4.9083E+03 

4.8959E+03 

20 

5.1382E+03 

5.0510E+03  4.9956E+03 

4.9547E+03  4.9223E+03 

4.8949E+03  4.8714E+03 

4.8509E+03 

4.8326E+03 

4.8163E+03 
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Fig.  5.5  Energy  efficiency  in  different  concentrations  of  input  solution 


5.2  Web  tension  control  system  simulation  results 


In  this  section,  simulations  of  the  web  tension  controller  will  be  provided. 
There  are  three  different  control  algorithms  to  be  compared  in  this  study.  A 
general  SMC,  and  two  SMCs  with  different  saturation  functions.  SMC_S  represents 
the  sliding  mode  controller  with  a  linear  saturation  function,  and  SMC_AS 
represents  it  with  an  automatically  adjusted  saturation  function.  Table  5.26  lists 
the  parameters  used  for  the  simulation,  and  table  5.27  is  the  control  gain  and  the 
control  goal  of  each  case.  Figs.  5.6,  5.9,  5.12,  5.15,  5.18,  5.21,  5.24  and  5.27  show 
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the  performance  of  the  sliding  mode  controller  in  different  cases.  As  indicated, 


there  are  some  issues  with  control  approach.  The  amplitude  of  chattering  of 
controlled  velocity  is  very  high,  and  the  tension  overshoot  is  about  15%  of  the 
desired  tension  setpoint.  From  figs. 5. 7,  5.10,  5.13,  5.16,  5.19,  5.22,  5.25  and  5.28 
show  the  performance  of  the  sliding  mode  controller  with  a  linear  saturation 
function  that  has  better  performance  than  the  general  SMC  especially  in  improving 
the  chattering  effect  in  controlled  velocity.  However,  the  tension  overshoot  still 
exists  (7%).  On  the  other  hand,  the  performance  of  SMC_AS  as  shown  in  Figs. 5. 8, 
5.11,  5.14,  5.17,  5.20,  5.23,  5.26  and  5.29  can  solve  the  chattering  and  overshoot 
problems.  In  addition,  the  transient  response  is  improved  and  the  settling  time  is 
about  250  times  less  than  the  other  two  methods.  For  the  SMC_S  controller, 
changes  in  the  web  tension  gain  can  improve  the  overshoot  problem,  but  the 
settling  time  will  still  be  significantly  greater  than  for  SMC_AS.  However  with 
SMC_S,  if  reduced  settling  time  is  necessary,  the  overshoot  problem  will  be  worse. 
Hence,  SMC_AS  is  a  better  solution. 

In  the  real  situation,  there  will  be  unexpected  disturbances  that  may  change 
the  real  tension.  Fig. 5. 30  and  Fig.  5.31  are  the  simulation  results  of  performance  of 
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the  SMC_S  and  SMC_AS  controllers  with  the  inclusion  of  30  dB  white  SNR  noise.  As 
illustrated,  both  of  these  methods  can  handle  this  situation  very  well.  Although 
there  is  some  chattering  caused  by  noise,  the  controllers  are  able  to  control  to  the 
required  set-point.  Comparing  the  performance  of  the  three  methods,  SMC_AS  is 
more  appropriate  for  this  study  for  the  following  reasons:  (1)  unexpected  high  web 
tension  can  cause  the  web  to  break,  and  (2)  even  if  the  material  does  not  fail, 
characteristics  of  web  material  may  change  from  this  deformation  and  this  can 
adversely  affect  the  fluid  dynamics.  Therefore,  for  our  study,  the  third  method  will 
be  better  choice. 
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Table  5.26  Pre-set  parameters  of  the  web  tension  control  system 


E 

Plastitic  of  elastic  modules 

396825(N/m^2) 

A 

Cross  section  Area  of  the  web 

0.0001(m^2) 

P 

Density  of  the  web 

1203.466(kg/m^3) 

Ljk 

Length  between  jth  and  kth  rollers 

3  (m) 

Lij 

Length  between  1th  and  jth  rollers 

1  (m) 

Lti 

Length  between  1th  roller  and  tail  pulley 

1.5  (m) 

w 

Width  of  the  web 

0.2  (m) 

Ri 

Radius  of  i-th  roller 

0.15  (m) 

Rj 

Radius  of  j-th  roller 

0.15  (m) 

Rk 

Radius  of  k-th  roller 

0.15  (m) 

Rm 

Radius  of  m-th  roller 

0.15  (m) 

Ri 

Radius  of  1-th  roller 

0.15  (m) 

Rt 

Radius  of  tail  pulley 

0.15  (m) 

h 

Thickness  of  the  web 

500  (pm) 

Ki 

velocity  constant  of  i-th  motor 

27.54  (N.m/volt) 

Kj 

velocity  constant  of  j-th  motor 

27.54  (N.m/volt) 

Kk 

velocity  constant  of  k-th  motor 

0.2391  (N.m/volt) 

Km 

velocity  constant  of  m-th  motor 

0.2391  (N.m/volt) 

Di 

Damper  coefficient  of  i-th  roller 

0.002(N.m.sec) 

Dj 

Damper  coefficient  of  j-th  roller 

0.002(N.m.sec) 

Dk 

Damper  coefficient  of  k-th  roller 

0.002(N.m.sec) 

Dm 

Damper  coefficient  of  m-th  roller 

0.002(N.m.sec) 

Table  5.27  Control  Parameters 


Case  I 

Case  II 

Case  III 

Case  IV 

Goals 

T=  IN, 

Vw  =1  in/min 

T=  IN, 

Vw=10 

in/min 

T=  IN, 

Vw  =20 

in/min 

T  =  ION, 

Vw  =1  in/min 

Tension  control  gain 

20 

20 

20 

20 

velocity  control  gain 

2 

2 

2 

2 

Case  V 

Case  VI 

Case  VII 

Case  VIII 

Goals 

T  =  ION, 

Vw=10 

in/min 

T  =  ION, 

Vw  =20 

in/min 

T  =  20N, 

Vw  =1  in/min 

T  =  20N, 

Vw  =20 

in/min 

Tension  control  gain 

20 

20 

20 

20 

velocity  control  gain 

2 

2 

2 

2 

Fig  .5.6  Performance  of  Sliding  mode  controller  (1  N  tension,  1  in/min 

transporting  velocity) 
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Fig.  5.7  Performance  of  SMC_S  (1  N  tension,  1  in/min  transporting  velocity) 
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Fig.  5.8  Performance  of  SMC_AS  (1  N  tension,  1  in/min  transporting  velocity) 
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Fig.  5.9  Performance  of  SMC  (1  N  tension,  10  in/min  transporting  velocity) 


108 


SICj 


tiiiie(sec) 


Fig.  5.10  Performance  of  SMC_S  (1  N  tension,  10  in/min  transporting 

velocity) 
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Fig.  5.11  Performance  of  SMC_AS  (1  N  tension,  10  in/min  transporting 

velocity) 
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Fig.  5.12  Performance  of  SMC  (1  N  tension,  20  in/min  transporting  velocity) 
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Fig.  5.13  Performance  of  SMC_S  (1  N  tension,  20  in/min  transporting 

velocity) 
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Fig.  5.14  Performance  of  SMC_AS  (1  N  tension,  20  in/min  transporting 

velocity) 
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Fig.  5.15  Performance  of  SMC  (10  N  tension,  1  in/min  transporting  velocity) 
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Fig.  5.16  Performance  of  SMC_S  (10  N  tension,  1  in/min  transporting 

velocity) 
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Fig.  5.17  Performance  of  SMC_AS  (10  N  tension,  1  in/min  transporting 

velocity) 
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Fig.  5.18  Performance  of  SMC  (10  N  tension,  10  in/min  transporting 

velocity) 
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Fig.  5.19  Performance  of  SMC_S  (10  N  tension,  10  in/min  transporting 

velocity) 
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Fig.  5.20  Performance  of  SMC_AS  (10  N  tension,  10  in/min  transporting 

velocity) 
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Fig.  5.21  Performance  of  SMC  (10  N  tension,  20  in/min  transporting 

velocity) 
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Fig.  5.22  Performance  of  SMC_S  (10  N  tension,  20  in/min  transporting 

velocity) 
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Fig.  5.23  Performance  of  SMC_AS  (10  N  tension,  20  in/min  transporting 

velocity) 
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tiffle(sec) 


Fig.  5.24  Performance  of  SMC  (20  N  tension,  1  in/min  transporting  velocity) 
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tiiiie(sec) 


Fig.  5.25  Performance  of  SMC_S  (20  N  tension,  1  in/min  transporting 

velocity) 
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Fg.  5.26  Performance  of  SMC_AS  (20  N  tension,  1  in/min  transporting 

velocity) 
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Fig.  5.27  Performance  of  SMC  (20  N  tension,  20  in/min  transporting 

velocity) 
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tiiiie(sec) 


tiiiie(sec) 


tiiiie(sec) 


Fig.  5.28  Performance  of  SMC_S  (20  N  tension,  20  in/min  transporting 

velocity) 
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Fig.  5.29  Performance  of  SMC_AS  (20  N  tension,  20  in/min  transporting 

velocity) 
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Fig.  5.30  Performance  of  SMC_S  controller  with  30  db  SNR  (20  N  tension,  20 

in/min  rewinding  velocity) 
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Fig.  5.31  Performance  of  SMC_AS  controller  with  30  db  SNR  (20  N  tension, 
20  in/min  rewinding  velocity) 
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5.3  Simulation  of  the  HDD  system 


As  discussed  previously,  the  main  goal  of  the  HDD  system  is  to  provide 
qualified  products  with  minimum  energy  for  the  production  of  high  quality 
fuel.  The  key  for  reducing  energy  loss  is  to  drain  as  much  water  as  possible 
before  the  wet  algae  enters  heater.  Hence,  the  fluid  dynamics  in  the  gravity 
well  and  capillary  zone  are  the  main  concerns  and  key  determining  factors.  In 
order  to  control  these  two  factors,  dominate  elements  such  as  input  solution 
concentration,  solution  temperature,  web  tension,  and  web  transport  velocity, 
are  required  to  be  adjusted  to  achieve  the  goal  of  saving  energy  and  providing 
as  much  dried  algae  as  possible  within  the  restrictions  of  the  material,  dryer 
system,  and  transport  system.  Fig.  5.32  is  the  diagram  for  the  system,  and 
table  5.28  provides  the  conditions  desired  for  the  simulation.  From  figs.  5.33~ 
5.36,  the  simulations  show  that  the  gravity  flow  rate  and  capillary  flow  rate 
are  very  stable,  and  the  energy  required  to  produce  one  pound  of  dried  algae 
when  there  exits  15  db  snr  white  noises  caused  by  web  tension  fluctuations. 
From  the  discussion  earlier,  the  performance  of  the  web  tension  controller  is 
proven  to  be  good,  and  it  is  proven  that  it  can  be  used  to  effectively  operate 
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the  HDD  system.  Even  as  unknown  effects  caused  by  tension  fluctuations  are 


introduced  to  the  system,  the  quality  of  the  yield  of  algae  can  still  be 
maintained.  To  sum  up,  energy  requirement,  and  yield  of  dried  algae  can  be 
maintained  at  desired  levels  even  when  low  or  high  white  Gaussian  noise 
disturbances  affect  the  process.  Also,  we  observe  that  the  performance  is 
good  for  different  web  velocity  and  web  tension  control  settings.  Hence,  we 
may  conclude  that  the  control  method  is  appropriate  for  the  HDD  system. 


Table  5.28  Simulation  cases 
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Fig.  5.32  System  control  module 
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Fig.  5.33  Simulation  of  HDD  system  (10  g/L,  10  N,  10  in/min,  15  db  SNR) 
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Fig.  5.34  Simulation  of  HDD  system  (10  g/L,  ION,  20  in/min,  15  db  SNR) 
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Fig.  5.35  Simulation  of  HDD  system  (10  g/L,  20N,  10  in/min,  15  db  SNR) 
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Fig.  5.36  Simulation  of  HDD  system  (20  g/L,  20N,  10  in/min,  15  db  SNR) 
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Chapter  6  Summary  and  Future  work 


In  this  thesis,  a  predictive  model  is  developed  for  HDD  system.  It  includes 
CFD  modeling,  roller-to-roller  web  tension  control  system,  and  heating 
control  system.  According  to  the  analysis,  energy  efficiency  and  dried  algae 
yield  are  highly  related  to  the  concentration  of  the  input  solution.  Modulating 
web  tension  can  be  used  to  reduce  the  energy  load  of  the  dryer.  If  conditions 
such  as  dryer  operating  settings,  web  velocity  limitations,  restrictions  on 
input  concentration,  and  production  cost  of  the  input  solution,  and  so  on,  can 
be  given,  the  desired  control  settings  can  be  determined  from  a  lookup  table. 
In  this  study,  a  general  dryer  is  chosen,  and  the  results  show  that  the 
operating  cost  for  dewatering  can  be  reduced  to  around  $1.5  dollars  per 
gallon  of  bio-diesel.  If  higher  concentrations  of  algae  can  be  used,  the 
dewatering  cost  can  be  further  reduced.  In  short,  the  method  provided  in  this 
study  offers  a  methodology  for  a  company  or  production  center  to  determine 
the  most  efficient  operation  of  the  harvesting  system  for  different  conditions. 
The  lookup  table  for  determining  operating  parameters  of  the  HDD  system 
can  reduce  the  trial  and  error  involved  with  the  control  of  the  system.  In 
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addition,  the  modeling  effort  is  an  effective  tool  for  the  design,  development. 


validation,  commercialization,  operation,  and  maintenance  of  the  HDD 
system. 

The  predictive  model  has  been  built,  and  it  seems  that  the  proposed 
HDD  system  can  be  operated  well  and  ride  through  some  unexpected 
disturbances,  although  the  control  system  has  not  been  tested  in  a  laboratory 
prototype.  Hence,  experiments  on  the  HDD  system  prototype  should  be  done 
in  the  future.  Besides,  the  verification  of  scaling  up  the  prototype  equipment 
is  required  after  a  fully  functional  control  system  has  been  implemented  and 
tested  in  the  HDD  prototype.  Moreover,  in  order  that  algae  can  be  a  market 
replacement  for  petroleum,  the  efficiency  of  the  production  process  has  been 
maximized.  Once  a  small-scale  prototype  is  built  and  the  HDD  is  ready  to 
transition  to  a  commercial  product  platform,  the  modeling  tools  developed  in 
this  work  will  also  transition  to  an  application  engineering  toolbox.  This  will 
serve  to  establish  optimal  system  designs  for  new  applications  and  to  reliably 
establish  the  expected  throughput,  energy  efficiency,  product  quality, 
operating  cost  and  economic  value.  The  core  elements  will  form  the  basis  for 
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an  adaptive  model-based  control  system  that  will  permit  the  system  to 


automatically  respond  to  system  changes  such  as  belt  wear,  algae  size 
variation,  biomass  dilution,  contaminants,  and  other  system  disturbances 
while  maintaining  reliable  and  efficient  production  levels.  In  addition, 
information  obtained  from  comparing  the  simulation  model  with  real-world 
system  performance  will  provide  an  effective  real-time  diagnostic  tool  that 
will  be  used  to  pinpoint  system  components  or  processes  that  are  degrading 
or  beginning  to  fail.  If  the  overall  production  cost  can  be  reduced  to  100 
USD/bbI,  algae-based  fuel  has  the  potential  to  be  a  competitive  alternative  to 
fossil  fuels. 
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